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ABSTRACT 
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The  purpose  of  this  research  was  to  Investigate  the  effect  of  tidal- 
period  Internal  wave  vertical  motion  on  oceanic  mixed  layer  dynamics, 
and  to  discern  the  effect  upon  the  diurnal  evolution  of  thermocllnes, 
as  observed  during  MILE  (mixed  layer  experiment).  Vertical  advectlon 
was  added  to  a  one-dimensional  bulk  model  of  the  mixed  layer  with  an 
assumed  linear  In  z  and  sinusoidal  In  time  dependence.  The  rate  of 
mixed  layer  deepening  was  therefore  due  to  the  combination  of  vertical 
motion  and  entrainment.  The  first  significant  result  was  the  finding 
that  the  Interaction  between  vertical  mixing  and  vertical  motion 
depended  upon  the  wave  frequency  and  Its  phase  relation  to  the  diurnal 
heating  cycle.  Second,  linear  and  non-linear  Interactions  of  the  wave 
Induced  vertical  motton  with  the  cyclical  boundary  conditions  can 
generate  two-dimensional  (x-z)  structure  In  the  near-surface  temperature 
field  of  an  Initially  horizontally  homogeneous  ocean  under  the  Influence 
of  horizontally  homogeneous  surface  boundary  conditions.  Finally,  this 
advectlve  Interaction  Increases  the  utility  of  the  mixed  layer  model  In 
single  station  forecasting. 
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I .  INTRODUCTION 


In  recent  years  much  study  of  Internal  waves  has  focused  on  the 
linear  and  non-linear  wave  Interactions.  Models  for  wave  prediction 
have  been  developed  and  tested  to  match  the  theory.  There  have  been 
attempts  to  measure  the  Internal  wave  fields  In  situ.  Some  of  the 
most  useful  of  these  have  been  the  observations  of  temperature  and 
pressure  field  variations  with  depth*  providing  a  basis  for  computa¬ 
tions  of  the  vertical  motions  and  Interactions  Involved.  These  studies 
and  analyses  of  data  demonstrate  the  Importance  of  the  Internal  wave 
field  In  upper  ocean  and  near  surface  dynamics. 

Roth  et  al.  (1981)  presented  clear  evidence  that  Internal  waves 
are  not  as  critical  In  an  analysis  of  the  dynamics  of  mixing  at  depth 
because  they  have  little  variability.  However,  near  the  surface, 
where  they  Interact  with  other  phenomena  such  as  cyclical  heat  flux, 
wind  stress,  and  other  boundary  conditions,  their  effect  can  be 
considerable  when  dealing  with  the  budgets  of  heat  and  momentum. 
Holloway  (1980)  pointed  out  that  Internal  waves  are  not  weak  waves 
(In  the  sense  of  weak  resonant  Interaction  theory)  and  have  a  high 
Impact  on  mixed  layer  dynamics  down  to  and  Including  turbulent  scales 
of  motion.  Plnkel  (1981)  showed  that  the  near-surface  wave  field  has 
a  profound  effect  on  the  vertical  velocity  field  In  the  upper  ocean. 

Some  mixed  layer  experiments  such  as  those  by  Linden  (1975)  and 
Kantha  et  al.  (1977)  have  shown  that  Internal  wave  radiation  may  play 
an  Important  role  In  controlling  the  rate  of  deepening  of  a  turbulent 


mixed  layer  In  a  stratified  fluid.  Displacements  Induced  by  the  turbu¬ 
lence  In  the  base  of  the  mixed  layer  served  to  form  Internal  waves 
which  removed  kinetic  energy  of  motion  from  the  mixed  layer.  Since 
this  energy  Is  also  Important  In  the  entrainment  process,  these  Induced 
Internal  waves,  of  relatively  high  frequency,  serve  to  retard  deepening 
of  a  mixed  layer.  Townsend  (1968)  has  derived  an  expression  for  the 
rate  of  generation  of  the  Internal  waves  and  their  energy  by  a  steady 
disturbance  field  which  Is  advected  relative  to  a  stratified  layer  at 
a  constant  speed.  The  theory  was  extended  by  Bell  (1978)  to  cover 
general  space-time  dependence  of  the  disturbance  field.  These  works 
Indicate  that  Internal  waves  could  be  a  means  of  retarding  the  growth 
of  mixed  layers,  through  the  removal  of  the  kinetic  energy  In  the  layers. 

Such  studies  all  Indicate  that  the  presence  of  Internal  waves  should 
be  taken  Into  account  as  a  possibly  Important  factor  In  the  mechanical 
energy  budget  which  controls  mixing.  It  Is  hypothesized  that  Internal 
waves  which  have  periods  comparable  to  the  dominant  time  scales  of  the 
atmospheric  forcing,  (surface  heat  flux  and  wind  stress)  will  be  the 
most  likely  to  result  In  an  Interaction  with  the  entrainment  mixing 
process  In  terms  of  a  positive  feedback  (Increased  mixing)  effect. 

A.  BACKGROUND  DATA  ANALYSIS 

During  the  period  August-September  1977,  a  mixed  layer  experiment 
(MILE)  was  conducted  at  Ocean  Station  P  (55°  N,  145°  W)  in  the  North 
Pacific.  The  depth  of  the  mixed  layer  varied  between  about  five  and 
thirty  meters  during  the  period  of  the  experiment.  Temperature  was 
continuously  observed  at  fixed  depths  for  a  period  of  three  weeks.  A 
spectrum  analysis  of  this  data  (See  Figure  10  In  later  section)  shows 
some  Important  features. 
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First,  a  major  spectral  peak  was  at  about  12.79  t  .02  hours,  with 
secondary  peaks  at  6.62  +  .13  and  3.94  ±  .06  hours.  Lower  order 
spectral  peaks  are  at  4.9  hours,  and  at  several  frequencies  having 
periods  under  3.5  hours.  The  peak  at  12.79  hours  Is  clearly  asso¬ 
ciated  with  the  semi-diurnal  tide  (12.45  hrs.),  with  the  other  two 
peaks  at  6.62  and  3.94  hours  probably  being  higher  harmonics.  At  all 
three  depths  of  the  analysis  It  Is  clear  that  the  semi-diurnal  internal 
tidal  wave,  with  its  harmonics,  Is  the  dominant  cause  of  vertical 
motion  at  the  base  of  the  mixed  layer. 

Second,  the  tidal  peaks.  Including  the  harmonics,  are  very  coherent 
vertically  throughout  the  depths  analyzeo.  They  are  present  In  all 
three  depths  analyzed,  below  the  layer,  and  In  the  mixed  layer,  with 
very  little  frequency  or  phase  shift. 

Finally,  the  spectral  peaks  at  other  than  tidal  frequencies  are 
not  vertically  coherent.  Also,  there  1s»a  large  variation  in  non-tldal 
spectral  densities  for  observations  above  and  below  the  thirty  meter 
depth— the  top  of  the  seasonal  thermocllne  at  this  time  of  year.  This 
suggests  that  something  has  happened  to  the  higher  frequency  waves 
near  the  surface  making  the  interpretation  of  their  effect  on  mixed 
layer  dynamics  hard  to  judge  from  simple  theories.  However,  at  the 
same  time  It  suggests  that  the  effect  of  tidal  internal  waves  may  be 
quite  easily  discerned. 

These  results  confirm  the  work  of  Plnkel  (1981)  whose  studies 
showed  that  of  the  wave  fields  In  the  upper  ocean,  the  tide  and  Its 
harmonics  have  the  greatest  vertical  coherence. 

Studies  by  Simpson  and  Paulson  (1979)  demonstrated  the  variability 
of  the  temperature  structure  In  the  upper  ocean  during  the  POLE  experiment. 


Shook  (1980)  showed  that  these  variations  were  a  response  to  vertical 
and  horizontal  advectlve  processes.  In  a  study  of  Pacific  Ocean  ship 
track  temperature  data*  Butler  (1981)  found  definite  cyclical  varia¬ 
tions  In  the  sea  surface  temperature  structure. 

From  these  studies  of  actual  data.  It  Is  obvious  that  (1)  the  tidal 
period  Internal  waves  and  their  harmonics  are  an  Important  source  of 
energy  In  the  upper  ocean,  and  (2)  surface  layer  temperature  patch¬ 
iness  exists  In  the  ocean,  and  may  In  part  be  generated  by  interaction 
of  mixing  and  Internal  wave  motion. 

B.  ADVECTION  IN  MODELLING 

To  date  little  has  been  published  concerning  the  interaction  of 
tidal-period  Internal  waves  and  entrainment  by  the  mixed  layer.  Be 
Szoeke  (1980)  has  discussed  advectlon  In  the  mixed  layer  as  a  result 
of  horizontal  variability  of  the  surface  wind  stress.  He  showed  that 
advectlon,  upwelllng  and  downwelllng,  was  an  Important  parameter  In 
mixed  layer  dynamics,  resulting  from  the  boundary  forcing  of  the  wind 
stress.  Cushman-Rolsln  (1981),  Garwood  et  al.  (1981),  and  Price  (1981) 
also  demonstrated  that  advectlon  can  play  an  Important  role  In  the 
dynamics  of  the  mixed  layer,  but  the  Interaction  with  Internal  gravity 
waves  Is  not  treated.  De  Szoeke  and  Rhlnes  (1981)  discussed  coastal 
upwelllng  at  length,  but  the  advectlve  motion  was  conceived  as  a 
necessary  factor  to  achieve  a  realistic  balance  for  surface  heating, 
transport,  and  mixing.  In  order  to  produce  observed  coastal  temperature 
fields.  Price  et  al,  (1978)  Included  vertical  advectlon  In  the  test 
of  a  mixed  layer  formulation  during  storm  conditions.  They  found  that 
the  phase  dependence  between  wind  stress  and  mixed  layer  deepening, 
and  the  mixed  layer  deepening  Itself  were  reasonably  simulated. 
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Within  this  discussion  of  a  United  amount  of  relevant  literature 
on  the  advectlve  process  It  has  been  shown  that  vertical  advectlon 
must  play  an  Important  role  In  storm  response*  coastal  circulation* 
and  fronts.  No  one  has  examined  the  role  of  Internal  wave  notion,  on 
the  tidal  scale*  In  the  open  ocean  where  the  geometry  and  boundary 
conditions  are  rather  one-dimensional. 

In  possibly  the  most  relevant  previous  effort,  Splgel  (1980) 
looked  at  the  coupling  of  Internal  wave  motion  and  entrainment  In  a 
two-layer  lake  model,  having  a  mixed  layer  which  Is  based  upon  the 
one-dimensional  turbulent  kinetic  energy  budget.  Time  scales  of 
internal  wave  motion,  wave  decay,  and  entrainment  were  widely 
separated.  Using  a  seiche  effect  In  a  basin,  he  demonstrated  that 
Internal  waves  did  In  fact  play  an  Important  role  In  the  entrainment 
zone. 

With  this  background  It  was  believed  that  the  time  has  came  for 
an  Investigation  of  the  effects  of  Internal  vertical  motion  In  the 
oceanic  mixed  layer.  The  motion  Is  there.  Observations  show  that 
It  dominates  Eulerlan  measurements  In  the  upper  thermocllne.  The 
Importance  of  this  motion  to  mixed  layer  evolution  needs  to  be 
established, 

C.  SCOPE  OF  THIS  STUDY 

The  scope  of  this  Initial  study  shall  be  confined  to  a  limited 
number  of  wave  frequencies.  There  is  Insufficient  time  to  concentrate 
on  all  the  types  of  Internal  vertical  motion,  and  all  of  the  possible 
couplings  and  Interactions  In  this  thesis.  The  purpose  of  this  Initial 
research  was  to  Investigate  the  effect  of  dominant  tidal-period  Internal 
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wave  vertical  Motion  on  oceanic  Mixed  layer  dynamics,  and  to  discern 
the  effect  upon  the  transient  thermocllne  as  observed  during  the  Mixed 
layer  experlMent.  The  effects  of  the  vertical  Motion  are  explored 
using  a  one-dimensional  bulk  model  of  the  Mixed  layer  originally  pro* 
posed  by  Garwood  (1976,1977)  and  expanded  here  to  Include  a  prescribed 
vertical  velocity. 

The  first  hypothesis  to  be  examined  Is  can  alteration  of  the 
existing  Mixed  layer  Model  be  reasonably  done?  Can  a  relatively 
simple  Mathematical  transformation  be  applied  to  Include  a  prescribed 
vertical  Motion  without  altering  the  basic  physics  of  the  entrainment 
process?  The  study  will  show  that  this  In  fact  can  be  done,  as  revealed 
by  a  careful  examination  of  the  Interactions  at  each  step  of  the  alter¬ 
ation  of  the  original  Model. 

A  second  hypothesis  to  be  examined  Is  that  vertical  notion  produced 
by  long  period  Internal  gravity  waves  for  which  the  Rossby  number  Is 
small,  through  Interaction  with  cyclical  boundary  conditions  (wind 
stress  and  heat  flux)  can  produce  two -dl Mens lonal  (x-z)  structures  In 
the  temperature  field  of  a  simple  Infinite  ocean  having  horizontally 
homogeneous  surface  boundary  conditions.  An  objective  will  be  to 
demonstrate  the  frequency  dependence  of  the  internal  vertical  motion 
upon  the  phases  of  the  cyclical  heat  fields,  and  the  wind  stress.  It 
will  be  shown  that  there  Is  a  non-linear  Interaction  between  super¬ 
imposed  diurnal  and  semi-diurnal  Internal  tidal  waves.  Finally,  this 
study  will  show  the  relative  strengths  of  the  resultant  two-dimensional 
structures,  and  relate  these  to  their  Importance  In  forming  an  improved 
thermodynamic  ocean  model. 
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II.  MIXED  LAYER  MODEL  AND  VERTICAL  MOTION  THEORY 
A.  MIXED  LAYER  MOOEL 

To  study  the  effects  of  Internal  vertical  motion  on  mixed  layer 
evolution,  the  Garwood  (1976*1977)  mixed  layer  model  was  used  In  a  one- 
dimensional  form.  Garwood  developed  an  ocean  mixed  layer  model  using 
the  Navler-Stokes  equation  of  motion  with  the  geostrophlc  component 
eliminated,  the  continuity  equation  for  Incompressible  water,  the  heat 
equation  from  the  first  law  of  thermodynamics,  the  conservation  of  salt 
equation,  and  a  linearized  equation  of  state.  The  entrainment  hypothe¬ 
sis  depends  upon  the  relative  distribution  of  turbulent  energy  between 
horizontal  and  vertical  components  and  Is  offered  as  a  plausible 
mechanism  for  governing  both  entrainment  and  layer  retreat. 

Garwood  suggested  that  planetary  rotation  Influences  the  dissipa¬ 
tion  of  turbulence  for  deeper  mixed  layers  and  enables  a  cyclical  steady 
state  to  exist  on  an  annual  basis.  Furthermore,  the  rate  of  entrainment 
for  the  stable  regime  cannot  be  a  simple  linear  extrapolation  of  the 
unstable  situations.  Unlike  the  atmospheric  boundary  case,  most  of  the 
solar  radiation  does  not  penetrate  the  layer.  Therefore,  downward 
turbulent  heat  flux  In  the  oceanic  boundary  layer  Is  as  Important  as 
the  upward  flux  during  the  course  of  both  diurnal  and  annual  cycles. 

The  non-linearity  of  the  Interface  entrainment  tendency  parameter, 
which  Is  greatest  for  stable  surface  boundary  conditions,  results  In 
a  modulation  of  the  long  term  trend  of  the  mixed  layer  depth  by  the 
diurnal  component  of  surface  heat  flux.  In  this  model,  buoyant  pro¬ 
duction  Is  somewhat  more  efficient  than  shear  production  as  a  source 
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of  energy  for  vertical  mixing  because  of  Its  unique  effect  on  the 
vertical  component  of  the  turbulent  velocity. 


The  governing  heat  equation  for  this  model  assumes  two  forms: 
me  for  the  mixed  layer,  the  other  for  below  the  mixed  layer.  For 
the  mixed  layer,  the  equation  has  the  form 


0) 


where  <  >  denotes  the  mixed  layer  average  and  Q  Is  the  net  heat  flux, 
a  product  of  downward  solar  radiation,  upward  turbulent  heat  flux  to 
the  surface,  and  surface  back  radiation.  Below  the  mixed  layer  the 
heat  equation  becomes 

dT  m  -d(wT)  ♦  aoQ  +  Kd*T 

i r  »  TcT  ^  (2) 

where  K  Is  the  eddy  viscosity,  and  «£(wT)/&z  Is  the  term  added  for  the 
vertical  motion. 

The  buoyancy  flux  at  the  entrainment  zone  Is  modelled  according  to 

-  -«77>“«y»  (3) 

which  Is  derived  from  the  local  turbulent  kinetic  energy  (TKE)  budget 
where 

z  ■  -h  (3a) 

In  the  entrainment  zone,  the  buoyant  damping  attributed  to  (wTTr)_h  Is 
assumed  to  be  balanced  by  the  convergence  of  TKE  flux.  This  convergence 
Is  modeled  as  a  function  of  the  TKE  kinetic  energy  components  and  the 
vertical  distance  over  which  this  energy  must  be  transported.  <wTwT> 
and  &  are  the  bulk  values  of  the  vertical  component  and  the  total  TKE, 

IB 


respectively.  The  prognostic  equations  for  these  mixed- layer  turbu¬ 
lence  variables  are  derived  using  the  bulk  second-order  closure  Methods 
of  Garwood  (1977).  and  are 

d  (h<E>)  -  "»«I  -  «qh(7TrLh  ♦  eghftTTLh  ♦  (7t\1 

-  (<£>**  fh)<i>  (♦) 


a  (h<wV>) 
at  2 


where 


•thftwT1)^  ♦  (7t\  1  ♦ 

2 

(<E>  -3<wrwS)<E>,/*  -  (<£>»»»fh)<g>  (5) 


u*  *  (r/^)1/2 

with  r  being  the  Magnitude  of  the  surface  wind  stress  and  p9  Is  the 
air  density,  assumed  to  be  a  constant,  and 

R*  «  OfhAT 

(Air+Av*) 


Is  the  bulk  Richardson  number.  The  change  In  the  total  TKE,  <T>,  Is  a 
function  of  wind  stress  shear  production,  entrainment  shear  production, 
bouyant  product Ion/ damping,  and  dissipation.  The  change  In  the  vertical 
component,  cwVV,  Is  a  function  of  bouyant  production/damping,  pressure 
redistribution,  and  dissipation.  The  quantities  AT,  Au,  and  Av  are  the 
jump  values  between  the  nixed  layer  (assumed  to  be  vertically  uniform), 
and  the  underlying  stable  waters.  These  bulk  TKE  equations  are  solved 
algebraically  assuming  quasi-steady  state  for  the  TKE  budgets 
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^L(h<T>)  -  o 

and 

d  (h<w*  *'>)  •  0 

at 

With  this  assuMptlon  the  computed  values  of  <w'w'>  and  <?>  are  In 
phase  with  the  surface  boundary  conditions— ■ the  wind  stress, feui,  and 
the  buoyancy  flux,  «g(w"T)Q. 

Once  the  entrainment  bouyancy  flux  Is  known  from  (3),  the  downward 
fluxes  of  momentum  associated  with  entrainment  at  the  base  of  the  mixed 


layer,  (wJu')-h  and  (wTvT)-h 

are  given  by 

- 

(6) 

-(w'v'Ljj  “ 

■,A*  -  -<w'T\.4s. 

(7) 

Here  the  entrainment  velocity  Is 

*e 

- 

• 

(3) 

and  the  entrainment  velocity  is  also  defined  as 

% 

"  If  + 

(9) 

If  w8  is  greater  than  zero.  If  wg  Is  not  greater  than  zero,  equations 
(10)  and  (11)  are  applied. 

1.  Boundary  Layer  Deepening  by  Entrainment 

The  entrainment  heat  flux  (viT)^  Is  determined  algebraically 
from  (3),  and  the  steady-state  forms  of  (4)  and  (5)  In  terms  of  given 
values  of  u*,  (?»TTr)0,  h  and  R^.  This  heat  flux  Is  then  applied  to 
the  given  temperature  profile.  The  new  h  and  new  momentum  and  temperature 
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profiles  are  unique  solutions  provided  that  (1)  the  mixed  layer  is 
homogeneous,  (11)  the  value  of  h  Is  Just  large  enough  to  prevent 
unstable  density  profiles  (dT£z  greater  than  zero),  and  (111)  heat 
Is  conserved. 

2.  Boundary  Layer  Shallowing 

Shallowing  occurs  when  there  is  adequate  vertical  turbulent 
energy  to  transport  heat  down  to  the  base  of  the  existing  mixed  Iqyer. 
The  depth  at  which  the  vertical  heat  flux  (w'T1)  vanishes  establishes 
the  new  layer  boundary,  and  will  be  the  base  of  the  new  mixed  layer. 

In  this  model,  this  occurs  when  (wTwT)  approaches  zero.  Then  (3)  Is 
no  longer  applicable,  and  the  steady  state  forms  of  (4)  and  (5) 
reduce  to 

O  «  *u?  -  (<£>*♦  fh)<i>  ♦  eqh(ffiA 

2  (10) 

n  -  «>»*  -  (<E>^*  fh)<E>  ♦  ogh(wrTr)A 

»  T  on 

Without  the  planetary  rotation  term,  or  If  the  downward  surface  buoyancy 
flux  Is  dominant,  h  Is  proportional  to  the  Obukhov  length  scale 

L  .  aj/«g(wT*)0 

If,  on  the  other  hand,  (wT )  •  0,  h  Is  proportional  to  u^/f.  In 

o 

general,  however,  the  depth  of  the  shallowing  mixed  layer  Is  a  function 
of  the  two  non-dimensional  parameters,  h/L  and  hf/u*. 

Heat  and  momentum  are  conserved  In  the  profile  specifications 


and,  because  the  mechanical  energy  budget  requires  conservation  of 
potential  energy  on  the  time  scale  of  the  turbulence,  a  numerical 
procedure  Is  used  that  also  preserves  potential  energy.  This  Is 


Intended  to  ensure  that  the  deepening  rate  Is  as  correct  as  possible 
when  the  mixed  layer  deepens  to  Its  former  position. 


B.  VERTICAL  MOTION  THEORY 


By  following  a  water  particle  as  It  Is  vertically  advected  by  a 
prescribed  velocity  field,  w(x,z,t),  the  distortion  of  any  temperature 
profile,  T(z),  after  a  small,  but  finite  length  of  time,  (At),  may  be 
computed  by  non-linear  theory. 

If  we  assume 


-  Az  coe(«t  +kx) 
0 


Is  prescribed,  then 

||  m  ooe(««t  +  kx) 

This  may  be  reformed  and  Integrated 


m  kx)dt 


(12) 


(13) 


where  a  water  particle  starting  at  time  at  position  z^  Is  advected 
to  position  z2  at  time  t2.  Or  following  any  constant  Isotherm 
Initially  at  position  zj,  to  z2  at  gives 


*,  •  t,exp{jL[»l«(«fc  “  eietwl,  +  kx)]} 

hO  I**) 

If  t  •  tj>  -  tj «  ,  then 


** 


*i  exp{A^i[ooeMt,  ♦  J)  ♦  kx)]} 


(15) 


which  Is  the  general  for*  approximated  by  the  progression  of  equations 
In  this  study  of  Internal  vertical  notion.  Internal  vertical  notion 
Is  applied  to  the  nixed  layer  model  through  equation  (10).  dh/dt* 
due  to  the  Internal  vertical  advectlon,  is  derived  from  equation 
(16). 

For  future  study*  the  salinity  budget  can  be  handled  In  a  like 
manner. 
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III.  EXPERIMENTS 


A.  PRODUCTS 

The  products  used  to  Identify  features  In  the  model  runs  and  used 
for  comparison  are: 

1.  Table  (Sample  Figure  la)  showing  the  hourly  change  In  the  mixed 
layer  primary  variables: 

H  Is  the  mixed  layer  depth  In  meters. 

T  Is  the  mixed  layer  temperature  In  degrees  centigrade. 

S  Is  the  mixed  layer  salinity  In  parts  per  thousand. 

TAU  Is  the  surface  wind  stress  In  dynes  per  square  centimeter. 

QNET  Is  the  net  heat  flux  out  of  the  mixed  layer  In  watts 
per  square  meter. 

WDD  Is  the  value  of  the  vertical  motion  parameter  affecting 
the  model  In  centimeters  per  second. 

DH  Is  the  change  In  mixed  layer  depth  due  to  the  vertical 
motion  for  that  time  period. 

2.  Contours  of  the  temperature  variation,  T(z,t),  (Sample  Figure  1) 
between  z  •  0  and  z  •  -60  meters  over  the  Initial  24  hours  of  the  model 
run.  Varying  resolution  Is  used  to  better  resolve  the  evaluation  of 
T(z,t).  Temperature  contours  are  In  degrees  centigrade. 

3.  Contours  of  temperature  variation  similar  to  2,  In  a  quasi  two- 
dimensional  mode.  This  is  accomplished  by  running  the  model  24  times, 
each  time  passing  the  prescribed  vertical  motion  through  the  24  hour 
cycle  of  the  model  at  a  different,  progressive  phase.  The  24  runs  are 
then  sorted  by  hour  giving  24  data  sets  of  24  Independent  points  with 

a  spatial  separation  of  kx  from  hour  one  through  hour  twenty-four. 
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It  Is  assumed  that  there  Is  a  two-dimensional  Internal  wave  In  x  and 
2  passing  through  this  ocean  system.  Heat  flux  and  wind  stress  are 
taken  to  be  the  same  at  all  points (  with  no  Interaction  between  any 
of  the  points  (advectlon  Is  limited  to  the  vertical). 

4.  Progressive  temperature  profiles,  T(z),  to  a  depth  of  60  meters 
(Sample  Figure  lc).  Only  the  first  profile  matches  the  horizontal 
temperature  scale  which  Is  degrees  centigrade.  The  remaining  profiles 
have  been  consecutively  offset.  The  accompanying  table  has  the  mixed 
layer  temperatures  In  degrees  centigrade,  and  the  mixed  layer  depth 

In  meters. 

5,  Progressive  temperature  profiles  as  In  4,  but  In  the  two- 
dimensional  mode  of  3.  Profiles  are  offset  from  the  first,  again,  to 
provide  an  insight  to  the  two-dimensional  nature  to  the  data. 

B.  STANDARDS 

Prior  to  the  examination  of  the  addition  of  vertical  motion  to  the 
mixed  layer  model  (MLM) ,  it  is  beneficial,  here,  to  examine  the  results 
of  running  the  mixed  layer  model  In  the  two  modes  which  will  be  used 
as  the  standards  for  comparison  with  the  vertical  motion  modelling 
experiments.  Both  standard  cases  are  runs  of  the  mixed  layer  model 
from  which  the  effects  of  evaporation  and  precipitation  have  been 
removed.  This  has  been  done  to  aid  In  the  interpretation  of  later 
results  by  limiting  the  number  of  Interaction  variables  to  a  manageable 
quantity.  The  first  standard  Is  a  run  of  the  MLM  in  the  primarily  wind- 
driven  deepening  mode.  The  second  standard  Is  for  the  MLM  In  a  cycli¬ 
cal  shallowing/deepening  response  to  a  strong  downward  heat  flux. 
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1.  Standard  1  (SI) 

This  run  was  Initialized  with  a  temperature  profile  having  a 
surface  and  mixed  layer  temperature  of  15.0°C.  At  5.0  ra  the  Initial 
temperature  decreases  linearly  with  the  form 

T(z)  .  T,  -  5.0  X  IO-4(Z-H)  (1?) 

Where  Z  Is  the  absolute  value  of  z  In  centimeters ,  and  H  Is  the  mixed 
layer  depth  In  centimeters.  This  profile  extends  to  200m.  The  result¬ 
ing  profile  Is  shown  In  Figure  2. 

The  surface  wind  stress  (TAU)  Is  a  constant  0.50  dynes  per 
square  centimeter  for  the  entire  run.  Net  heating  reaches  a  maximum 
of  302.5  watts  per  square  meter  at  hour  six  (local  noon).  Net  heat 
loss  after  sunset  Is  a  constant  100.8  watts  per  square  meter. 

With  no  vertical  motion  the  layer  deepens  continuously  over 
the  two-day  period  of  the  model  run  (Figure  1),  as  expected  for  pri¬ 
marily  wind-driven  mixing.  In  examining  the  first  twenty-four  hour 
period  It  Is  seen  that  the  mixed  layer  deepens  from  an  Initial  depth 
of  4.5  m  to  a  depth  of  14.7  m.  The  layer  cools  from  15°C  to  14.767*C 
over  the  period,  with  a  short  period  of  mixed  layer  temperature 
Increase  of  six  hours  related  to  that  time  when  the  net  heat  flux 
Is  sufficient  to  overcome  the  cooling  attributed  to  the  entrainment 
of  cooler  water  into  the  layer. 

The  salinity  In  the  layer  remains  constant  because  there  is 
no  precipitation  or  evaporation  allowed  which  could  change  the  surface 
salinity,  and  the  below  layer  salinity  Is  assumed  to  be  a  constant 
34.00  part  per  thousand.  Therefore,  there  are  no  sources  of  water 
with  a  salinity  other  than  34.00  parts  per  thousand. 


27 


Depth  u» 


34 


Figure  2.  Initial  temperature  and  salinity  profiles  for  Case  SI. 
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The  downward  solar  radiation,  Q$,  Increases  from  sunrise  (hour  1) 
to  noon  (hour  6),  and  then  decreases  to  zero  at  sunset  (hour  12.)  A 
constant  Is  assumed  for  the  net  upward  turbulent  heat  flux  to  the  atmos¬ 
phere,  plus  back  radiation  of  100.8  watts  per  square  meter  for  the 
entire  day. 

The  one-day  response  to  this  cycle  can  be  seen  In  the  contours 
(Figure  1b).  The  mixed  layer  deepens  continuously  over  the  first  24 
hours.  Below  the  layer  temperature,  at  any  given  depth  not  penetrated 
as  the  layer  deepens ,  remains  constant  over  time.  In  the  layer  the 
temperature  Is  seen  to  cool  slightly  Initially,  then  warm  to  a  maximum 
value  at  about  hour  nine  (1500  local  time),  and  then  cool  again  for 
the  remainder  of  the  period. 

In  the  T(z)  graph  this  sequence  Is  seen  again  (Figure  lc). 

Here  it  Is  seen  that  all  of  the  slopes  of  the  profiles  below  the  layer 
are  the  same  from  hour  one  to  twenty-four.  Indicating  that  with  no 
vertical  motion,  the  characteristics  below  the  layer  will  not  change 
with  time.  The  Included  table  shows  the  relationship  of  the  temper¬ 
ature,  In  degrees  centigrade,  and  the  layer  depth.  In  meters,  and 
Indicates  the  rate  and  magnitude  of  the  deepening. 

The  rate  of  deepening  (Table  1)  reflects  the  surface  boundary 
conditions  at  any  given  time  as  well  as  layer  depth  and  stratification 
below.  As  the  run  begins,  the  layer  deepens  quickly  in  response  to 
the  surface  wind  stress  and  the  short  distance  over  which  it  must  act, 
as  seen  In  equations  (3)  and  (4).  At  hour  one  the  layer  deepens 
beginning  at  a  shallow  value  where  the  wind  flux  has  a  high  Impact 
from  equation  (4),  during  a  period  where  the  heat  flux  Into  the  layer 
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Is  a  minimum.  As  the  heat  flux  Increases  and  the  layer  deepens  the 
wind  stress  Is  relatively  less  effective  In  deepening,  and  the  entrain¬ 
ment  now  effected  by  equation  (S)  slows  until  hour  nine  when  the  solar 
radiation  goes  to  zero,  and  the  deepening  rate  Increases.  This  Increase 
Is  relatively  small  because  of  the  relatively  large  temperature  jump 
(AT)  In  the  entrainment  zone  by  this  time.  After  hour  eleven  (1800 
local  time)  the  minor  heat  loss,  small  wind  stress  and  increased  depth 
Involved  all  act  to  slowly  decrease  the  deepening  rate. 


Table  1.  Deepening  Rate  In  Meters  Per  Hour  for  Case  SI 


Hour 

Hour 

Deepening 

Rate 

1 

1.7613 

13 

0.3296 

2 

1.0227 

14 

0.3376 

3 

0.7271 

15 

0.3287 

4 

0.4877 

16 

0.3017 

5 

0.4242 

T7 

0.3083 

6 

0.3221 

18 

0.3115 

7 

0.2941 

19 

0.2759 

8 

0.2964 

20 

0.2815 

9 

0.2904 

21 

0.2872 

10 

0.2917 

22 

0.2703 

11 

0.3291 

23 

0.2560 

12 

0.3580 

24 

0.2607 

2.  Standard  2  (S2) 

This  run  Is  begun  with  the  same  Initial  layer  temperature  as 
for  standard  1.  Here,  however,  the  Initial  mixed  layer  depth  Is  taken 
to  be  49,5  meters.  At  50.0  meters,  the  profile  again  assumes  the 
linear  variability  of  equation  (18),  extending  to  200.0  meters 
(Figure  3). 

All  other  parameters  are  the  same  as  those  of  SI  except  for 
the  prescribed  surface  heat  flux.  The  net  heat  flux  has  increased  to 
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Figure  3.  Initial  t— parature  and  salinity  profiles  for  Case  S2. 

3\ 


a  maxi nun  of  504.2  watts  per  square  meter  by  assuming  an  extremely 
strong  solar  heating  maximum  (Q$)  of  600  watts  per  square  meter  at 
local  noon. 

In  this  case  the  layer  depth  remains  at  49.5  meters  for  the 
first  hour  (Figure  21a)  because  the  wind  stress  Is  not  strong  enough, 
coupled  with  the  available  heat  energy,  to  cause  further  entrainment 
at  this  depth,  and  the  solar  heating  Is  not  yet  strong  enough  to  cause 
(wTwr)  to  go  to  zero  so  that  equations  (10)  and  (11)  are  applied.  By 
the  second  hour  the  net  downward  surface  heat  flux  and  the  surface 
wind  stress  are  sufficient  to  cause  (w1 w1")  to  go  to  zero,  and  a  new  h 
Is  formulated  through  the  solution  of  equations  (10)  and  (11),  forming 
a  new  surface  mixed  layer  with  a  depth  of  about  18  m.  The  layer  then 
continues  to  shallow  until  hour  6  (local  noon)  in  response  to  the 
Increased  downward  heat  flux.  Stirring  by  the  wind  stress  Is  relatively 
small  as  Its  effect  Is  overcome  by  the  buoyancy  flux  associated  with 
net  surface  heat  flux.  The  layer  temperature  Increases  as  long  as 
there  Is  heating  sufficient  to  overcome  the  minor  cooling  due  to  the 
entrainment  of  cooler  below  layer  water  which  resumes  at  local  noon, 
but  Is  not  strong  until  hour  10  (1600  local  time).  At  that  time  the 
layer  cools  due  to  the  combined  effects  of  the  surface  heat  loss,  and 
the  entrainment  of  the  cooler  below  layer  waters. 

The  salinity,  again,  remains  constant  because  there  Is  no 
source  of  water  of  a  salinity  other  than  34.00  parts  per  thousand. 

As  shown  by  the  contours  (Figure  lib)  the  original  layer  of  50  meters 
shallows  to  hour  6,  not  forming  a  distinct  layer,  since  equations 
(10)  and  (11)  are  applied  at  each  hourly  step,  until  hour  7.  The 
transient  thermocllne  prior  to  hour  7  Is  weak.  This  Is  seen  more 
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clearly,  perhaps,  on  the  T(z)  profiles  (Figure  11c).  The  originally 
sharp  Interface  at  50  meters  becomes  diffuse  with  time  due  to  small 
(Itj.  *  0.2  cm^/s)  below-layer  diffusion.  This  Is  clearly  seen  In  the 
T(z)  profiles,  and  It  Is  the  reason  for  the  non-horizontal  contours 
near  the  50  m  depth.  The  variability  of  the  temperature  structure 
below  the  newly  formed  surface  mixed  layer  Is  a  function  of  the  shal¬ 
lowing  adjustment  by  equations  (10)  and  (11)  as  applied  from  hour  two 
through  six. 

Again  It  Is  evident  In  the  T(z)  profiles  that  the  profiles 
keep  the  same  slope  below  the  layer  except  for  the  originally  sharp 
Interface  at  50  meters  where  diffusion  of  the  layer  occurs  with  time. 
This  diffusion  area  Is  also  seen  to  hold  at  a  relatively  constant 
depth  with  time. 

The  T(z,t)  evolution  reflects  the  turbulent  kinetic  energy 
budget  for  this  problem  (Table  2).  The  Initial  decrease  In  h  from  a 
deep  layer  to  a  shallow  layer  Is  clearly  the  result  of  the  surface 
heating  and  the  surface  wind  stress  acting  to  form  a  new,  warmer, 
shallower  surface  layer.  Energy  Is  Increased  In  the  surface  mixed 
layer  due  to  excessive  solar  insolation  during  the  day.  The  Obukhov 
length  1$  reduced  from  hour  two  through  six,  and  the  newly  forming 
mixed  layer  adjusts  accordingly.  By  hour  seven  the  solar  Insolation 
begins  to  reduce.  As  this  happens,  the  Obukhov  length  increases, 
allowing  the  newly  formed  layer  to  deepen.  Deepening  does  not  keep 
up  with  the  Obukhov  length  because  of  the  damping  due  to  the  energy 
transport  which  must  take  place  during  the  mixing  process.  By  hour 
11  the  net  heat  flux  Is  negligible  decreasing  to  a  constant  surface 
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loss  after  sunset.  The  relative  magnitude  of  the  depth  Increase 
(dh/dt)  with  time  when  compared  with  that  of  SI  Is  due  to  the  strength 
of  the  temperature  jump  (AT)  at  the  layer,  and  the  available  energy 
In  the  layer.  In  SI  a  AT  of  the  order  of  .4°C  was  encountered, 
while  In  S2,  due  to  the  slow  formation  of  the  new  surface  layer,  a 
AT  of  the  order  .08°C  Is  the  maximum  encountered.  This  allows  for 
more  rapid  deepening.  The  rate  of  deepening  Increases  after  hour  16 
(2200  local  time)  because,  as  seen  In  the  T(z)  profiles  (Figure  11c) 
the  AT  at  the  entrainment  zone  Is  decreasing  with  time  due  to  the  net 
upward  surface  heat  flux. 


Table  2.  Deepening  Rate  In  Meters  Per  Hour  for  Case  S2. 


Hour 

mmm 

Deepening 

Rate 

1 

0.0152 

13 

0.8464 

2 

-31.4079 

14 

0.8317 

3 

-1.9247 

15 

0.8252 

4 

-3.1646 

16 

0.8235 

5 

-0.9575 

17 

0.8278 

6 

-0.2935 

18 

0.8404 

7 

0.3766 

19 

0.8536 

8 

0.3445 

20 

0.8879 

9 

0.4944 

21 

0.9186 

10 

0.5383 

22 

0.9516 

11 

0.6793 

23 

0.9870 

12 

0.8371 

24 

1.0256 

C.  EXPERIMENTS 

The  addition  here  to  the  existing  theory  of  Garwood  (1976-77)  Is 
approached  In  a  step-by-step  procedure,  gradually  Increasing  the  com¬ 
plexity  of  the  problem.  A  quick  reference  for  this  progression  can 
be  seen  In  Table  3,  A  basic  assumption  that  was  made  In  all  of  the 
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Table  3.  Initial  conditions*  boundary  conditions,  and  progression 
of  vertical  notion  Inputs  for  all  experiments. 


Par  all  Caaaa: 

1.  Mac  surface  haat  laaa  ta  carat aat  at  100. t  watts/a2. 

2.  tarfaca  wind  atsaaa  la  eoaataat  at  O.S  dyaes/ea2. 

0.  The  laitlal  layer  fpaaatraa  la  1S.0*C. 

4.  Tka  aal laity  la  eoaataat  ta  200  a  at  S4.00*/oo. 

5.  Mo  fraelpltatloa  or  evaporatloa  la  alloaad  Marine  the  rua. 

4.  The  vortical  aotloa  la  oapoaantlally  decreasing  to  the  surface.  aotioa 

la  hall  at  34.0  ca/hr  at  10  a,  escape  for  Caao  IVb,  where  it  la  10.0  ca/hr  at 

10  a. 


Caao 

I.  C. 

1. 

C. 

Vertical  Notion 

(wetta/h*) 

Profile 

MJ> 

<■) 

SI 

400.00 

Fif.  2 

4.S 

Naae 

S3 

000.00 

FI*.  3 

40. S 

-More 

la 

SI 

01 

01 

Coaataat  dawawslllng  during  aodel  rua. 

lb 

a 

02 

02 

Ie 

U 

01 

01 

01 

Ccaetaat  upselling  during  aodel  rua. 

» 

02 

02 

I  la 

01 

01 

01 

Fttlaad  dovnwelling,  12  hour  cycle.  Phaee 
shifted  to  peak  at  hours  5,  IS,  and  24. 

lib 

S3 

02 

02 

lie 

01 

01 

01 

Ftolaed  upeelliag,  12  hour  cycle.  Phase 
shifted  to  peak  at  hours  5,  IS,  and  24. 

ru 

02 

02 

02 

Ilia 

01 

01 

01 

Slagle  diurnal  wave  passage.  Phase  shifted 
to  give  earl—  downward  aotioa  at  all  hours. 

Illb 

02 

02 

02 

Hears  1,  7,  IS,  and  19  discussed. 

m 

02 

02 

02 

CeatlaOoas  diurnal  wave  systsa  passage.  Saae 

phase  shifts  as  Casa  III.  Hours  6,  12,  is. 

24  dilCBINd* 

m 

S3 

02 

02 

Coatlaeows  seal -diurnal  wave  systsa  passage. 

Haas  shifted  as  la  Case  III.  Runs  with  peek 

■ettaas  at  hours  6  sad  1(,  sad  12  sad  24  are 

Vh 

02 

02 

03 

Sapertapoeed  diurnal  sad  seal -diurnal  continuous 

wavs  systsa  passage,  loth  waves  la  phase. 

Input  phase  shifted  as  la  Case  III. 

Vk 

02 

02 

02 

As  la  Casa  Va,  but  phases  of  diurnal  sad  seal- 

aavaa  varied  relative  to  oae  aaether. 
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-sues. 


work  was  that  the  vertical  notion  moved  the  water  up  and  down,  but  did 
not  directly  Induce  turbulence  (l.e.,  mixing).  This  Is  a  valid 
assumption  If  no  shear  Is  associated  with  this  ration. 

The  cases  of  progressively  complicated  vertical  ration  examined 
were: 

1.  Linearly  (In  z)  varying  vertical  ration,  constant  over  the 
entire  period. 

a.  Upwelllng 

b.  Downwelllng 

2.  Linearly  varying  vertical  ration  pulsed  In  time, 
a.  Short  dally  upwelllng/ downwelllng  events. 

3.  Linear  Internal  wave  motion. 

a.  Passage  of  a  single  wave  pulse  through  the  experimental 
point  over  the  24  hour  dally  cycle. 

4.  Non-linear  wave  variations  coupled  with  diurnal  and  semi¬ 
diurnal  boundary  conditions  (heating,  wind). 

5.  Non-linear  response  tosuper  Imposed  waves  of  different  periods. 
Diurnal  and  semi-diurnal  tidal  waves  are  combined  and  phase- 
shifted  to  determine  non-linear  Interactions. 

1.  Case  I 

Both  upwelllng  and  downwelllng  are  realistically  Induced  with 
the  realization  that  In  a  mixed  homogeneous  layer  the  vertical  dis¬ 
placement  due  to  a  linear  w(z)  can  be  shown  to  be  exponential,  going 
to  zero  at  the  surface.  Mathematically  we  know  that 

-  WM  (18) 
if 

w(z)  - 

where  D  Is  a  scaling  depth.  It  follows  that 
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or 


dz 

TT 


ii 

z 


Az 


-  H>- 


Integration  fron  tj  to  tg  yields 

*a  .  z, 


(19) 


which  Is  the  equation  upon  which  the  following  work  Is  based,  and  which 
agrees  with  equation  (16)  derived  earlier. 

The  temperature  and  salinity  profiles  must  reflect  "correct" 
variation  If  equations  (16)  and  (19)  are  correct.  Upwelllng  and/or 
downwelllng  will  not  cause  Isollnes  to  break  the  surface  or  each  other. 
Mixing  must  be  done  to  Induce  this.  Therefore,  above  the  mixed  layer 
depth  Tngw  and  must  equal  TQld  and  SQld  when  there  Is  vertical 
motion  without  mixing  or  surface  exchange.  Below  the  entrainment  depth 
the  profiles  must  be  altered  to  reflect  the  type  of  vertical  notion 
(upward  or  downward)  In  accordance  with  equation  (16).  Figure  4  shows 
the  expected  result  from  upwelllng  on  a  temperature  profile.  The  mixed 
layer  depth  has  decreased  fron  time  one  to  time  two  (after  upwelllng), 
and  the  slope  of  the  profile  below  the  mixed  layer  has  been  decreased. 

In  all  the  cases  the  AT  at  the  entrainment  zone  has  been  preserved, 
a.  Constant  Downwelllng  (Case  la)  (Figure  12) 

The  Initial  conditions  used  are  the  same  as  those  for  SI 
(See  Table  3).  A  constant  downward  vertical  motion  Is  Introduced  giving 
the  vertical  motion  parameter,  WDD,  (HDD  ■  A  t/D)  a  value  of  .03600, 
which  equates  to  a  value  of  vertical  motion  of  .01  cm/s  or  36  cm/hr  at 
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Temperature  profile  variation  from  tj  to  tg  under  upwelllng  conditions. 


ten  meters.  This  Is  a  rather  large  value,  but  It  demonstrates  the  Inter¬ 
actions  well.  The  vertical  motion  Is  applied  beginning  the  first  hour, 
and  Is  not  applied  as  the  Initial  condition  on  the  temperature  and 
salinity  profiles. 

The  result  of  this  downward  motion  Is.  as  expected,  a  more 
rapid  deepening  of  the  surface  mixed  layer  (Figure  12a)  when  compared 
to  SI.  Comparing  the  depth  change  due  to  the  mixing  in  this  case  (Table 
4)  to  the  Si  mixing,  the  same  sequence  of  events  as  described  for  SI  Is 
seen,  and  the  values  are  comparable.  Indicating  that  the  physics  of  the 
model  has  not  been  altered  by  the  addition  of  the  vertical  motion  para¬ 
meter  In  the  downwelllng  mode.  It  is  noted  that  the  effect  of  the 
vertical  motion  quickly  becomes  greater  than  the  effect  of  mixing  due 
to  the  large  value  of  vertical  mixing  assigned  In  this  run  of  the  model. 


Table  4.  Deepening  Rate  in  Meters  Per  Hour  for  the 
Mixing  and  Vertical  Motion  of  Case  la. 


Hour 

Mixing 

Deepening 

Rate 

Vertical  Motion 
Deepening 

Rate 

Hour 

Mixing 

Qeepenlng 

Rate 

Vertlcle  Motion 
Deepening 
Rate 

1 

1.7613 

0.2295 

13 

0.3201 

0.5729 

2 

0.9705 

0.2735 

14 

0.3169 

0.6055 

3 

0.6606 

0.3077 

15 

0.3134 

0.6392 

4 

0.4759 

0.3365 

16 

0.3102 

0.6740 

5 

0.3535 

0.3618 

17 

0.3068 

0.7099 

6 

0.2945 

0.3858 

18 

0.3034 

0.7471 

7 

0.2444 

0.4089 

19 

0.3017 

0.7855 

8 

0.2282 

0.4323 

20 

0.2983 

0.8252 

9 

0.2394 

0.4569 

21 

0.2949 

0.8663 

10 

0.2513 

0.4828 

22 

0.2914 

0.9087 

11 

0.2793 

0.5108 

23 

0.2881 

0.9526 

12 

0.3219 

0.5413 

24 

0.2094 

0.9952 

The  sequence  Is  vividly  displayed  In  the  contours  (Figure  12b) 
where  It  Is  seen  that  the  temperature  contours  deepen  as  time  progresses 
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as  one  would  expect  In  a  downwelllng  case.  The  surface  heating  due 
to  the  Incoming  solar  radiation  Is  noted  again,  followed  by  cooling 
as  the  sun  sets  and  heat  is  then  lost  at  the  surface  and  colder  water 
is  entrained  at  the  bottom  of  the  mixed  layer.  It  is  also  noted  that 
In  only  twenty-four  hours,  the  vertical  motion  added  here  has  deepened 
the  layer  to  almost  twice  the  value  of  SI. 

The  T(z)  plot  (Figure  12c)  demonstrates  the  deepening 
Involved  and  the  infects  on  the  temperature  profile  below  the  mixed 
layer.  Note  the  Increased  slope  (dT/dz)  with  time,  as  wanner  near¬ 
surface  temperatures  (water)  are  carried  downward  by  the  downwelling. 
Also,  one  can  see  the  effect  of  the  preservation  of  the  AT  profile  at 
the  entrainment  zone  from  the  mixing  step  through  the  vertical  motion 
progression  In  the  model  where  at  times  there  Is  an  apparent  discon¬ 
tinuity  In  the  profile.  This  does  not  appear  to  cause  a  problem  In 
the  execution  of  the  model,  and  will  be  seen  In  many  of  the  following 
T(z)  plots.  One  can  further  see  that  the  vertical  motion  addition  has 
caused  a  slight  decrease  In  the  jump  strength  of  AT  at  the  layer  over 
that  In  SI.  This  can  be  attributed  to  the  reduced  mixing  for  the 
larger  mixed  layer  depths.  Where,  with  no  downwelllng  a  AT  of  1 .0°C 
may  have  occurred  In  18  meters,  It  now  occurs  at  a  greater  depth. 

This  effect  Is  transmitted  Into  the  entrainment  zone  as  this  below 
layer  water  Is  mixed  In,  causing  the  AT  to  be  reduced, 
b.  Constant  Downwelllng  (Case  lb)  (Figure  13) 

The  Initial  conditions  are  those  of  S2  (See  Table  3).  A 
constant  downward  motion  Is  Introduced  equal  to  that  used  In  the 
previous  case. 
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The  result  of  the  downward  motion  Is  again,  as  expected. 


a  more  rapid  deepening  of  the  layer  after  the  new  surface  mixed  layer 
Is  fonaed  through  the  surface  wind  stress  and  heat  flux  as  applied  by 
equations  (10)  and  (11)  (Figure  13a).  Comparing  the  depth  change  due 
to  the  mixing  here  (Table  5)  to  the  S2  mixing  shows  that  the  values 
are  of  the  same  order,  which  shows  that  the  energy  relations  of  the 
original  model  are  not  Impaired  by  the  addition  of  vertical  motion 
In  this  mode  of  the  model  either.  Here  It  is  noted  that  mixing  drives 
the  production  of  a  new  surface  mixed  layer  beginning  at  hour  2  (0800 
local)  (Figure  13b).  The  shallowing  due  to  the  surface  heating  Is 
dominant  until  hour  seven  (1300  local)  when  the  downwelllng  becomes 
the  dominant  factor  until  hour  nine.  From  hour  nine  to  hour  fifteen 
the  effects  of  mixing  and  downwelllng  are  of  the  same  order.  After  hour 


Table  5.  Deepening  rate  In  meters  per  hour  for  the 
Mixing  and  Vertical  Motion  In  Case  lb. 


Hour 

Mixing 

Deepening 

Rate 

Vertical  Motion 
Deepening 

Rate 

Hour 

Mixing 

Deepening 

Rate 

Vertical  Motion 
Deepening 
Rate 

1 

0.0152 

1.8150 

13 

0.8093 

0.6986 

2 

-33.4052 

0.6571 

14 

0.7343 

0.7511 

3 

-2.5282 

0.5885 

15 

0.7445 

0.8059 

4 

-3.7245 

0.4735 

16 

0.7131 

0.8616 

5 

-1.4432 

0.4392 

17 

0.7029 

0.9190 

6 

-  .7286 

0.4286 

18 

0.6768 

0.9775 

7 

0.1233 

0.4489 

19 

0.6716 

1.0379 

8 

0.3275 

0.4773 

20 

0.6628 

1.1002 

9 

0.4266 

0.5104 

21 

0.6420 

1.1641 

10 

0.5787 

0.5504 

22 

0.6316 

1.2299 

11 

0.6446 

0.5942 

23 

0.6224 

1.2978 

12 

0.8003 

0.6986 

24 

0.6120 

1.3678 

fifteen  the  downwelllng  again,  slowly,  becomes  the  dominant  factor  for 
layer  deepening. 
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It  Is  seen  that  the  layer  temperature,  again.  Immediately 
begins  to  warn  (Figure  13c),  and  In  fact,  while  the  dominant  factor  Is 
the  large  amount  of  surface  heating,  the  layer  warms  at  exactly  the 
same  rate  as  In  S2.  After  hour  ten  the  layer  again  begins  to  cool,  but 
because  of  the  deepening  of  the  mixed  layer  due  to  downwelllng,  the 
mixing  of  the  surface  wind  stress  Is  reduced.  This,  along  with  the 
fact  that  a  deeper  layer  has  available  more  mass  of  water  to  dilute 
the  effect  of  the  entrained  water,  results  In  less  cooling  over  the 
remaining  twelve  hours. 

Because  of  the  downwelllng  effect  on  diffusing  the  temper¬ 
ature  gradient  It  Is  also  seen  (Figure  13b)  that  the  entrainment  zone 
(layer  boundary)  remains  rather  diffuse  until  hour  eight,  while  In  S2 
It  was  sharpened  by  hour  six.  This  helps  explain  the  magnitude  of  the 
deepening  effect  in  the  period  from  hour  nine  to  fifteen,  though  the 
depth  of  the  entrainment  zone  is  increasing  faster  than  that  In  S2  and 
reaches  a  difference  of  five  meters  from  that  In  S2  by  hour  fifteen. 

The  AT  across  the  bottom  of  the  layer  for  the  downwelllng  case  is 
slightly  less  than  that  In  S2,  allowing  the  mixing  to  be  a  more  impor¬ 
tant  factor  for  deepening  in  this  case.  The  following  decrease  In  the 
mixing  Is  then  due  to  the  ever  Increasing  depth  through  which  the  mix¬ 
ing  must  take  place.  The  wind  stress  must  act  over  distances  that  are 
twice  as  large  as  those  In  S2,  and  although  the  AT  at  the  layer  depth 
Is  not  as  great,  the  effect  of  the  depth  Is  quite  notice#*^. 

The  original  deep  layer  Is  seen  to  begin  to  diffuse,  but 
Is  quickly  carried  below  the  60  meter  level  of  the  graphs.  As  stated 
above,  the  new  surface  mixed  layer  begins  to  show  a  significant  AT 
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at  hour  eight.  Extremely  close  examination  of  the  T(z)  profiles  for  S2 
and  the  profiles  for  this  run  does  show  a  smaller  AT  at  the  entrainment 
zone  for  this  case  during  the  period  from  hour  seven  to  hour  sixteen. 

After  hour  16  there  Is  no  dlscernable  difference  In  the  AT  In  the  profiles, 
c.  Constant  Upwelllng  (Case  Ic)  (Figure  14) 

The  Initial  conditions  are  those  of  SI.  A  constant  upward 
vertical  motion  Is  Introduced  beginning  at  hour  one  giving  the  vertical 
motion  parameter,  MOD,  a  value  of  .03600,  or  about  36  cm/hr  at  ten  meters. 
Again,  this  Is  high,  but  It  Is  done  with  the  realization  that  later  work 
will  be  with  wave  motion  where  this  will  not  be  an  unrealistic  value, 
and  using  It  here  provides  a  good  benchmark  for  the  effects  produced  by 
this  motion. 


Table  5.  Deepening  Rate  In  Meters  Per  Hour  for  the 
Mixing  and  Vertical  motion  of  Case  Ic. 


Hour 

'Mixing 

Oeepenlng 

Rate 

Vertical  Motion 
Oeepenlng 

Rate 

Hour 

Nixing 

Deepening 

Rate 

Vertical  Motion 
Deepening 
Rate 

1 

1.7613 

-.2214 

13 

0.4239 

-.3026 

2 

1.1031 

-.2526 

14 

0.4528 

-.3079 

3 

0.6746 

-.2675 

15 

0.4845 

-.3142 

4 

0.6103 

-.2796 

16 

0.4774 

-.3199 

5 

0.4954 

-.2873 

17 

0.2909 

-.3189 

6 

0.3417 

-.2892 

18 

0.2871 

-.3178 

7 

0.3184 

-.2902 

19 

0.3046 

-.3173 

8 

0.3097 

-.2909 

20 

0.3223 

-.3175 

9 

0.3133 

-.2917 

21 

0.3407 

-.3183 

10 

0.3254 

-.2929 

22 

0.2596 

-.3163 

11 

0.3435 

-.2947 

23 

0.2739 

-.3148 

12 

0.3938 

-.2982 

24 

0.2891 

-.3138 

The  result  of  the  addition  of  this  upward  motion  Is,  as 
anticipated,  a  slowing  of  the  deepening  rate  of  the  mixed  layer  (Figure 
14a)  when  compared  to  SI.  In  fact,  due  to  the  size  of  the  vertical 
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notion  parameter,  an  effect  Is  produced  which  a 1 Most  holds  the  mixed 
layer  at  a  constant  depth.  This  Is  displayed  more  dramatically  after 
two  model  days  have  been  executed,  and  can  be  seen  to  reach  a  near 
equilibrium  condition  after  four  to  five  days.  Comparing  the  deepening 
rates  due  to  the  effect  of  the  mixing  In  the  model  and  the  vertical 
motion  in  this  case  (Table  5)  to  the  SI  mixing  It  can  be  seen  quite 
clearly  that  the  effects  described  for  the  mixed  layer  model  standard, 
SI,  are  still  present.  From  equations  (1)  and  (3)  the  slightly 
Increased  magnitude  of  the  cooling  of  the  mixed  layer  by  entrainment 
and  the  subsequent  heating  are  attributed  to  the  reduced  depth  over 
which  the  wind  stress  and  heat  flux  must  be  distributed.  Again,  It 
can  be  seen,  that  the  model  TKE  balance,  equations  (3-5),  Is  not 
violated  by  the  addition  of  this  upwelllng  motion. 

The  layer  temperature  Initially  falls  (Figure  14a),  then 
rises  qntll  1400  local  time.  It  then  falls  again  In  response  to  the 
mixing  and  the  surface  heating  cycle,  as  In  SI.  The  heating  of  the 
layer  is  not  as  great  here  as  In  the  downwelllng  case  (Figure  13), 
because  the  entrainment  Is  sufficiently  enlarged  due  to  the  smaller  h, 
and  because  cooler  water  becomes  available  for  entrainment  more  quickly 
due  to  the  upwelllng.  This  Is  more  clearly  seen  In  comparing  the 
contoured  fields,  T(z,t),  (Figure  14b)  of  this  case  with  those  of  the 
downwelllng  case  (Figure  13b).  By  the  end  of  the  twenty-four  hour 
period  this  has  resulted  In  a  very  strong  temperature  jump, AT,  at 
the  entrainment  zone,  as  Is  seen  best  In  the  T(z)  profiles  (Figure  14c). 
It  Is  seen  that  AT  grows  with  time. 

The  strengthening  of  the  jump,  and  the  lifting  of  the 
cooler  water  Is  demonstrated  clearly  by  the  contours,  and  It  Is 
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clearly  evident  the  contours  not  only  lift,  but  also  merge  In  time. 
Indicating  the  AT  grows  with  time.  This  effect  Is  seen  In  the  T(z) 
profiles  by  the  Increased  gradient  of  the  profiles  with  time. 

Equilibrium  Is  reached  In  this  case  as  the  jump  gets  too 
large  for  the  mixing  to  have  an  effect  greater  than  the  effect  of 
the  vertical  motion.  Equation  (8)  shows  that  dh/dt  vanishes  as  we 
approaches  w_^. 

d.  Constant  Upwelllng  (Case  Id)  (Figure  15) 

The  Initial  conditions  applied  are  Identical  to  those  of 
$2.  A  constant  upward  vertical  motion  is  Introduced  at  hour  one  of 
about  36  cm/hr  at  ten  meters,  as  before. 

The  net  result  of  the  vertical  motion  after  24  hours  Is  a 
reduced  layer  depth  In  comparison  to  the  case  having  no  vertical  motion 
(S2).  The  total  reduction  In  depth  of  six  meters,  however.  Is  not  as 
great  as  the  previous  case  (Ic).  After  two  days  It  Is  seen  that  th» 
mixed  layer  Is  approaching  an  equilibrium  depth  of  near  eleven  meters 
(Figure  15a). 

The  deepening  rates  (Table  6),  when  compared  to  those  of 
S2,  demonstrate  that  the  basic  model  Is  still  preserved  In  this  run. 
From  hours  two  through  five  the  model  Is  In  the  shallowing  mode,  and 
the  new  layer  depth  Is  established  by  equations  (10)  and  (11),  while 
there  Is  a  massive  Influx  of  solar  radiation  to  the  sea  surface.  The 
Increased  deepening  rate  due  to  the  mixing  seen  from  hour  twelve 
through  twenty-three  Is  attributed  to  the  effect  of  equations  (1), 

(4)  and  (5).  The  mixed  layer  Is  maintained  at  a  shallow  level  by  the 
upward  vertical  motion,  therefore,  the  mixing  takes  place  over  a  much 
smaller  distance  than  in  S2. 
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Table  6.  Deepening  Rate  In  Meters  Per  Hour  for  the 
Nixing  and  Vertical  Motion  of  Case  Id. 


Hour 

'  '  TffxTng 
Deepening 
Rate 

vertical  Motion 
Deepening 

Rate 

Hour 

Mixing 

Deepening 

Rate 

Vertical  Motion 
Deepening 
Rate 

1 

0.0152 

-1.7509 

13 

0.9465 

-.4751 

2 

-29.4675 

-0.6470 

14 

0.9470 

-.4917 

3 

-1.3404 

-0.5767 

15 

1.0159 

-.5103 

4 

-2.6147 

-0.4639 

16 

1.0611 

-.5298 

5 

-0.5166 

-0.4292 

17 

1.1449 

-.5515 

6 

0.1312 

-0.4186 

18 

1.2301 

-.5755 

7 

0.3085 

-0.4148 

19 

1.3043 

-.6013 

8 

0.4578 

-0.4163 

20 

1.4087 

-.6298 

9 

0.5885 

-0.4224 

21 

1.5577 

-.6626 

10 

0.6775 

-0.4314 

22 

1.6638 

-.6980 

11 

0.7113 

-0.4413 

23 

1.3743 

-.7219 

12 

0.9078 

•0.4578 

24 

0.6909 

-.7208 

Further,  there  Is  a  concentration  of  energy  In  this  shallower 
layer.  The  great  variance  In  the  last  hour  Is  a  direct  result  of  the 
magnitude  of  the  upwelllng,  which,  by  hour  four,  has  set  up  a  large 
gradient  in  the  temperature  (Figure  15b).  By  hour  twenty-four  a  gradient 
has  been  set  up  with  such  a  large  slope  that  mixing  Is  becoming  very 
difficult.  Continued  upward  motion  of  this  magnitude,  without  shallow¬ 
ing  as  Is  seen  on  the  second  day,  would  quickly  make  It  extremely 
difficult  to  entrain  much  water  Into  the  mixed  layer. 

The  progression  In  the  change  of  the  mixed  layer  temperature 
(Figure  15a)  shows  that  as  the  model  Is  shallowing,  governed  by  equation 
(10)  and  (11),  the  same  adjustment  takes  place  as  is  seen  In  S2.  The 
temperatures  produced  are  exactly  those  seen  In  S2.  This  Is  another 
Indication  that  the  model  Is  performing  as  before.  Since  we  Is  zero, 
the  upwelllng  has  no  effect  on  the  setup  of  the  new  layer.  If  the  H 
values  had  been  listed  lanedlately  after  the  mixing  phase  of  the  compu¬ 
tation,  this  could  be  seen  more  clearly.  After  the  layer  again  begins 


to  mix  downward,  however,  the  upwelllng  begins  to  have  some  effect  on 
T  and  h  due  to  the  upward  transport  of  the  entrainment  zone.  Though 
the  layer  Is  one  to  two  meters  shallower  than  In  S2,  the  temperature 
of  the  layer  Is  the  same  as  that  In  the  deeper  standard  to  hour  ten. 

At  hour  ten  the  heat  flux  of  equation  (1)  has  been  reduced  sufficiently 
so  that  wind  mixing  again  becomes  the  dominant  deepening  force.  At 
this  point  the  upwelllng  has  moved  water  cooler  than  that  In  S2  at  the 
same  time  to  the  entrainment  zone.  The  change  In  entrainment  rate  due 
to  the  upwelllng,  the  entrainment  of  cooler  waters,  and  the  net  surface 
heat  loss  after  hour  eleven  combine  to  produce  a  cooling  of  .147*C  from 
the  maximum  warm  temperature,  as  compared  to  .109*C  for  S2.  For  the 
mass  of  water  Involved,  and  the  relatively  short  time  Involved,  this 
Is  a  significant  difference. 

The  effect  of  the  upward  vertical  motion  Is  perhaps  best 
shown  In  the  T(z)  diagram  (Figure  15c),  where  It  Is  seen  that  the  model 
Initially  adjusts,  through  equations  (10)  and  (11),  to  a  new  surface 
layer  which  slowly  strengthens,  in  terms  of  the  jump,  AT,  as  It  en¬ 
trains.  As  the  process  continues  the  profiles  below  the  layer  rise 
by  upwelllng  and  the  Initial  deep  Interface  is  smoothed  by  diffusion 
until  It  disappears.  By  hour  twenty-four  a  sufficiently  stable  pro¬ 
file  has  been  set  up  to  make  deepening  of  the  layer  more  difficult. 

This  Indeed  is  the  case  In  the  second  day  of  the  run  (Figure  15a)  as 
the  layer  reaches  a  temporary  equilibrium  It  Is  temporary*  In  that 
once  the  profile  has  been  reached  which  makes  deepening  sufficiently 
difficult  through  equations  (4)  and  (5),  the  dominant  force  becomes 
the  upwelllng,  w_h  of  equation  (9). 
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The  contours  (Figure  15b)  further  demonstrate  the  effect 
of  the  upward  lifting  of  the  cooler  water,  as  the  contours  tend  to 
converge  with  time.  The  Initial  Interface  Is  quickly  lost,  as  the 
nodel  smooths  this  Into  Its  uplifted  profile. 

2.  Case  II 

The  next  step  of  this  progressive  series  of  tests  Is  to  explore 
the  effects  of  a  simple  pulsed  motion  varying  In  time.  This  Is  accom¬ 
plished  by  varying  the  Input  motion  In  a  sinusoidal  pattern.  By  using 
the  equation 

»0(t)  m  AA  sin  (2*1,  +  kx) 

**  (20) 

as  the  particular  form  of  equation  (12)  applied  for  this  case,  where 
AA  is  the  velocity  scale  factor  such  that  AA  »  Az/D  is  scaled  for  an 
appropriate  motion  at  100  meters,  and  limiting  the  value  of  w0(t)  to 
only  values  In  one  direction  (+  or  -) ,  and  zero,  a  pattern  of  motion 
can  be  produced  such  as  that  seen  In  Figure  5.  A  pulse  of  vertical 
motion  of  twelve  hour  duration  Is  produced,  followed  by  a  twelve  hour 
period  of  no  motion.  To  determine  the  relationship  of  this  motion 
with  the  cyclical  boundary  conditions  one  need  only  shift  the  phase 
of  the  pulse  relation.  This  will  cause  the  motion  to  peak  at  varying 
selected  times  during  any  twenty-four  hour  period  run. 

a.  Pulsed  Downwelllng  (Case  Ila)  (Figure  16,  17  and  18) 

The  Initial  conditions  are  the  same  as  those  of  SI.  Three 
different  phase  relationships  of  the  downwelllng  to  the  cyclical 
boundary  conditions  are  presented.  In  the  first  (Figure  16)  the 
maximum  of  the  downwelllng  occurs  180*  out  of  phase  with  the  heating 


48 


S'' 


cycle  (that  period  of  the  model  defined  as  night,  or  the  cooling  cycle). 
In  the  second  (Figure  17),  the  downwelllng  occurs  during  the  last  seven 
hours  of  night,  and  during  the  first  five  hours  of  daylight.  The  max* 
Imun  downwelllng  occurs  six  hours  before  the  heating  maximum.  Finally, 
In  the  third  case  (Figure  18),  the  downwelllng  occurs  In  phase  with 
the  maximum  heating  cycle.  In  all  three  cases,  the  downwelllng  event 
Is  a  twelve  hour  pulse  with  a  maximum  value  of  36  cm/hr  at  ten  meters. 

(1)  180*  Out  of  Phase  (Figure  16).  For  this  case  It  is 
seen  that  the  model  evolves  exactly  as  does  SI  until  hour  thirteen 
when  downwelllng  begins  (Figure  16a).  All  the  fields  produced  are 
exactly  the  same  as  In  SI  until  hour  thirteen  for  which  the  evolution 
has  been  discussed. 

After  hour  thirteen  the  downwelllng  begins  to  trans¬ 
port  the  surface  water  downward,  and  causes  a  decrease  In  the  jump,  AT, 
strength  of  a  few  hundredths  of  a  degree.  This  combined  effect  causes 
the  deepening  of  the  mixed  layer  to  increase  here  over  that  of  SI  by 
about  four  meters  over  the  twenty-four  hour  period.  Additionally,  it 
decreases  the  change  of  temperature  of  the  mixed  layer  because  the  layer 
Is  Initially  deeper,  thereby,  reducing  mixing,  and  Increasing  the  dilu¬ 
tion  of  the  entrained  waters,  and  by  causing  entrainment  of  water  which 
Is  not  as  cool  as  that  of  SI,  over  the  same  period.  These  effects  are 
apparent  In  the  contours  (Figure  16b),  and  In  the  T(z)  profile  (Figure 
16c).  The  downwelllng  effect  Is  seen  most  dramatically  In  the  contours, 
where  the  downward  pulse  Is  evident,  along  with  the  described  temper¬ 
ature  variations. 

(2)  Maximum  Downwelllng  Before  Heating  (Figure  17).  Here, 
the  downwelllng  Is  applied  during  the  last  six  hours  of  the  mixing. 
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and  during  the  first  five  hours  of  the  heating  cycle.  The  six  hours 
of  downwelllng  (hiring  the  last  six  hours  of  night  serve  only  to  deepen 
the  layer  In  the  manner  previously  described.  The  first  five  hours  of 
doMMelllng  coupled  with  the  heating  and  wind  stress  through  equation 
(9),  produces  a  noted  change  In  the  temperature  fields  (Figure  17c) 
from  that  of  SI.  This  coupling  leads  to  an  increase  of  the  layer 
temperature  at  the  peak  heating  period  of  .008°C  for  this  run  over 
that  of  the  previous  (8T  Increase),  and  over  that  of  the  standard,  SI. 
After  twenty- four  hours  this  difference  has  decreased  to  .005°C  as 
compared  to  the  previous  run,  and  by  two  days  the  difference  In  temp¬ 
eratures  Is  negligible,  though  the  layer  depth  is  slightly  shallower 
after  two  days  for  this  run.  By  only  partially  coupling  the  effects 
of  downwelllng  and  heating,  a  shallower  mixed  layer  has  been  produced 
with  somewhat  warmer  temperatures. 

The  surface  temperature  evolution  and  the  below  layer 
evolution  are  best  revealed  by  the  contoured  fields  (Figure  17b).  It 
Is  apparent  here  that  another  phenomenon  Is  present.  As  the  downwell¬ 
lng  coammnces  the  near  layer  contours  do  In  fact  diverge  as  expected. 
The  deeper  contours,  however.  Initially  narrow  as  the  pulse  forms, 
then  widen  to  the  below  layer  temperature  spreading  effect  described 
previously.  This  convergence  during  the  peak  motion  period  Is  a  direct 
result  of  the  mathematical  treatment  of  the  below  layer  temperature 
profiles  through  equation  (16). 

(3)  Coupled  Heating  and  Downwelllng  (Figure  18) .  Here 
the  downwelllng  peaks  one  hour  before  the  maximum  surface  heat  flux. 

The  combination  of  downwelllng  and  wind  stress,  equation  (9)  through 
equations  (16}  and  (8),  are  sufficiently  strong  to  cause  the  effect  of 
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the  heat  flux,  equation  (1),  to  be  negligible,  and  the  layer  deepens 
to  almost  three  meters  more  than  that  of  SI  by  hour  eleven.  Also, 
the  layer  temperature  has  been  Increased  significantly  over  SI  due 
to  the  downwelllng  effect  on  the  layer  which  causes  the  entrainment 
of  less  water  than  previously  during  the  period  of  maximum  surface 
heating.  By  combining  the  downwelllng  with  the  heating,  a  maximum 
temperature  change  of  .2U#C  for  SI  has  been  reduced  to  .162#C  for 
this  run  over  the  first  twenty-four  hours.  This  has  also  produced  a 
mixed  layer  which  Is  .06*C  warmer  and  three  meters  deeper  than  SI 
after  only  twenty-four  hours.  This  reflects  a  substantial  change  In 
the  energetics  of  the  system,  and  may  be  a  significant  mechanism  for 
the  generation  of  horizontal  temperature  structure.  After  the  pulse, 
the  layer  Is  mixed  by  the  action  of  the  wind  stress  In  the  presence 
of  a  smaller  AT— reduced  by  the  downwelllng  accompanied  by  a  small 
downward  surface  heat  flux.  This  explains  why  the  deepening  Is  not 
as  great  after  the  downwelllng  pulse  as  would  be  expected  If  one  con¬ 
sidered  only  the  reduced  jump  strength. 

In  combination  these  three  phased  runs  demonstrate 
the  effect  of  the  Interaction  of  the  downwelllng  and  the  diurnal 
heating  cycle.  In  all  cases  a  relative  heating  of  the  mixed  layer 
Is  seen  because  of  the  depth  Increase.  The  deeper  layer  reduces  the 
effect  of  the  mixing,  and  reduces  the  entrainment.  As  the  downwelllng 
event  couples  with  the  heating,  the  final  mixed  layer  depth  Is  reduced, 
over  that  of  the  uncoupled  cases,  and  greater  surface  temperatures 
(layer  temperatures)  are  achieved.  The  rudiments  of  a  pattern  In 
surface  features  which  could  be  Induced  over  a  large  ocean  by  slajriy 
varying  the  coupling  of  these  two  basic  effects  Is  now  seen. 
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b.  Pulsed  Downwelllng  (Case  lib)  (Figures  19,  20  and  21) 

The  Initial  conditions  are  those  of  S2.  The  magnitudes 
of  the  downwelllng,  and  the  phase  relations  are  the  same  as  the 
previous  case. 

(1)  180*  Out  of  Phase  with  Heating  (Figure  19).  Because 
the  downwelllng  event  occurs  after  the  surface  heating  has  stopped  (at 
night),  the  model  has  gone  through  the  shallowing  mode,  equations  (10) 
and  (11),  and  produced  a  new  surface  mixed  layer  exactly  as  In  S2  dur¬ 
ing  the  daylight  period.  The  effect  of  the  downwelllng  Is  seen  at  hour 
fourteen  when  mixed  layer  deepening  Is  first  observed.  The  cumulative 
effect  of  the  downwelllng  Is  to  deepen  the  layer  by  about  five  meters 
more  than  In  $2  In  the  twenty-four  hour  period  (Figure  19a).  The  layer 
temperature,  due  to  the  downwelllng  effect,  remains  .011°C  wanner  after 
the  period.  Destruction  of  the  original  layer  at  50  m  Is  caused  by 
diffusion,  and  the  model's  mathematical  smoothing  in  the  production  of 
new  temperature  contours  due  to  the  downwelllng  In  accordance  with 
equation  (16). 

(2)  Maximum  Downwelllng  Before  Heating  (Figure  20).  The 
shallowing  of  the  model  through  equations  (10)  and  (11)  during  the 
surface  heating  period  overrides  the  downwelllng,  so  that  by  hour  five 
this  run  Is  essentially  the  same  as  S2,  with  only  a  slight  temperature 
difference.  Comparing  the  contours  (Figure  11b  and  Figure  20b)  shows 
that  this  run  Is  the  sane  as  S2  until  the  downwelllng  pulse  begins  at 
hour  nineteen,  except  the  early  downwelllng  has  deepened  the  layer 
and  moved  the  cooler  water  away  from  the  entrainment  zone,  relative 

to  the  temperatures  In  S2.  The  layer  then  deepens  due  to  wind  stress 
mixing,  and  downwelllng.  The  net  result  of  this  coupled  action  Is  a 
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layer  about  four  meters  deeper  than  S2  after  twenty-four  hours*  but 
with  a  temperature  essentially  equal  to  that  of  S2. 


(3)  Coupled  Heating  and  Oownwelllng  (Figure  21).  Here 
the  shallowing  mode  is  still  the  dominant  factor,  with  the  downwelllng 
playing  only  a  minor  role  In  that  It  deepens  the  newly  forming  layer, 
and  moves  the  relatively  cooler  water  away  from  the  entrainment  zone. 
Because  the  model  Is  In  the  shallowing  mode,  as  we  approaches  zero, 
the  new  value  of  h  Is  dependent  on  equations  (10)  and  (11).  as  Is  the 
temperature  of  the  new  mixed  layer  through  equation  (1).  Since  the 
downwelllng  adjustment  Is  done  after  the  evaluation  of  the  shallowing 
h,  the  downwelllng  does  deepen  the  newly  formed  layer,  but  no  temper¬ 
ature  adjustment  Is  made.  This  results  In  a  layer  that  appears  to  be 
deeper  than  In  S2,  and  yet  maintains  the  same  temperature.  The  net 
result  of  these  interactions  after  twenty-four  hours  Is  a  layer  that 
Is  one  meter  deeper  than  In  S2,  and  ,012°C  wanner.  Here  the  coupling 
effect  was  negligible. 

In  combination  these  three  runs  of  the  model  In  the 
downwelllng  mode  with  S2  Initial  conditions  showed  the  opposite  pro¬ 
gression  of  those  runs  using  the  SI  Initial  conditions.  The  maximum 
variations  of  temperature  and  mixed  layer  depth  occurs  In  the  case  of 
180*  out  of  phase  downwelllng  with  the  heating  cycle,  where  a  deeper, 
wanner  mixed  layer  Is  formed  during  the  twenty-four  hour  period.  In 
the  case  of  the  SI  Initial  conditions  the  maximum  deviations  occurred 
during  the  coupled  run. 

C,  Pulsed  Upwelllng  (Case  Sic)  (Figure  22,  23  and  24) 

The  Initial  conditions  are  those  of  Si.  The  phase 


relationships  used  are  the  same  as  In  the  pulsed  downwelllng  case,  as 
are  the  magnitudes  of  the  effects. 


(1)  180*  Out  of  Phase  with  Heating  (Figure  221.  The 
effect  here  Is  i.jgllglble  until  hour  thirteen,  and  this  run  is  exactly 
the  same  as  SI  until  then.  At  hour  thirteen  the  upwelllng  begins  to 
increase  the  magnitude  of  AT.  Colder  water  Is  moved  up  toward  the 
entrainment  zone  below  the  layer,  and  the  layer  Itself  is  moved  upward, 
as  the  wind  stress  acts  to  deepen  the  layer  by  mixing.  At  hour  twenty- 
one  It  Is  seen  that  the  combination  of  the  upwelllng  and  the  magnitude 
of  AT  are  great  enough  to  actually  shallow  the  layer  for  that  hour. 

The  net  effect  of  the  upwelllng  in  this  case  Is  to 
create  a  mixed  layer  that  Is  about  2.6  meters  shallower  than  in  SI 
over  the  twenty-four  hour  period,  and  .03°C  cooler  (Figure  22a).  Also, 
there  Is  a  maximum  8T  In  the  mixed  layer  of  .241°C,  where  mixing  only 
produced  a  8T  of  .211°C  In  SI.  This  has  the  potential  of  producing 
an  observable  surface  temperature  effect. 

The  effect  of  contour  convergence  during  the  pulse 
mentioned  In  the  downwelllng  case  Is  seen  here  again  (Figure  22b). 

The  overall  effect  of  the  upwelllng,  to  converge  the  contours  from 
before  to  after  upwelllng.  Is  seen  as  expected. 

(2)  Maximum  Upwelllng  Before  Heating  (Figure  23).  In 
this  case  the  upwelllng  Immediately  begins  to  lift  the  entrainment 
zone  against  the  effect  of  the  wind  stress,  causing  an  Increase  In 
the  effect  of  the  mixing,  and  reducing  the  water  available  for 
dilution  of  the  entrained,  cold  water.  In  equation  (9),  w_h  is  of 
opposite  sign  to  that  of  wft,  reducing,  or  changing  the  sign  of  dh/dt. 
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This  results  In  cooler  layer  temperatures  over  the  same  period  as 
described  for  SI ,  and  shallower  mixed  layer  depths.  Numerically, 
the  net  result  of  tills  effect  over  the  twenty-four  hour  period  Is 
to  produce  a  layer  that  Is  2.S  m  shallower  than  that  of  SI,  and  ,Q4*C 
cooler  (Figure  23c).  The  maximum  ST  of  the  layer  over  the  period  Is 
,240*C,  or  about  the  same  as  In  the  previous  upwelllng  case. 

(3)  Coupled  Heating  and  Upwelllng  (Figure  24).  The 
coupling  of  the  upwelllng  and  the  heating  cycle  produce  a  mixed  layer 
that  Is  only  1.9  meters  shallower  than  that  of  SI  In  the  twenty-four 
hour  period,  but  a  layer  that  is  much  cooler.  The  reduction  of  the 
layer  depth  during  the  heating  period  causes  the  maximum  temperature 
to  be  .013*C  colder,  and  causes  the  layer  temperature  at  the  end  of 
the  twenty-four  hour  period  to  be  .Q71*C  colder  than  SI,  a  significant 
change.  The  maximum  ST  of  the  layer  for  this  period  Is  .269°C. 

As  In  the  downwelllng  case  with  the  shallow  initial 
layer  conditions,  the  maximum  surface  variations  In  temperature  Is  In 
the  coupled  upwelllng  case. 

d.  Pulsed  Upwelllng  (Case  I Id)  (Figure  25  and  26) 

The  Initial  conditions  are  those  of  S2.  The  phase  rela¬ 
tionships  used  are  the  same  as  those  In  the  previous  three  model  runs. 
Magnitudes  of  the  upwelllng  are  the  same  as  in. the  previous  run. 

(1)  180*  Out  of  Phase  with  Heating  (Figure  25).  The 
shallowing  proceeds  as  In  S2  for  the  first  six  hours,  and  a  new  layer 
Is  formed  (Figure  25b).  At  hour  thirteen  the  upwelllng  begins  and 
several  effects  are  noted.  First,  the  deepening  lamed lately  slows, 
relative  to  that  In  S2,  as  the  upwelllng  lifts  the  layer,  and 
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Increases  the  AT  at  the  entrainment  zone.  Second,  the  layer  temperature 
decreases  faster  relative  to  that  In  S2,  due  to  the  shallowing  of  the 
layer,  and  entrainment  of  relatively  cooler  water  Into  the  layer. 

Finally,  there  Is  an  effect  not  previously  observed.  The  Initial  layer, 
which  Is  being  diffused  In  the  previous  runs  is  now  reinforced  during 
the  upwelllng  (Figure  25c),  and  a  two-layer  system  Is  formed  over  the 
first  twenty-four  hours. 

The  net  result  of  this  case  Is  a  surface  layer  that  Is 
2.4  meters  shallower  and  .022°C  cooler  than  S2  over  the  period.  This 
leads  to  a  slightly  Increased  surface  signature  over  the  non-upwelllng 
case.  Because  of  the  reinforcement  of  the  Initial  layer  there  is  also 
a  potential  for  a  substantially  recognizable  T(z)  signature  for  this  case. 

(2)  Maximum  Upwelllng  Before  Heating  (Figure  26).  The  shal¬ 
lowing  and  upwelllng  combine  to  make  this  a  relatively  uninteresting 
case.  Though  the  large  upwelllng  has  some  initial  effect,  the  fact 
that  we  goes  to  zero  forcing  the  model  to  be  driven  by  equations  (10) 
and  (11)  soon  overcomes  the  upwelllng,  and  by  hour  four  this  case  is 
Identical  to  S2  and  the  previous  run  with  the  exception  that  the  up¬ 
welllng  has  advected  cooler  water  closer  to  the  entrainment  zone. 

The  difference  In  water  temperature  and  mixed  layer  depth  between 
this  case  and  S2  are  negligible  by  hour  eighteen.  At  hour  nineteen 
upwelllng  begins  again.  Increasing  the  relative  effect  of  the  wind 
stress  mixing,  causing  Increased  entrainment  of  cool  water.  The  net 
effect  of  the  upwelllng  and  the  wind  stress  mixing  over  the  depth 
involved  Is  to  produce  a  layer  that  Is  2.2  meters  shallower  than  in 
S2,  but  at  only  .004*C  cooler  after  twenty-four  hours.  This  would 


not  produce  a  noticeable  effect  at  the  surface.  The  two- layer  type 
system  seen  In  the  previous  case  could*  again*  provide  a  signature 
to  this  sort  of  upwelllng  event. 

(3)  Coupled  Heating  and  (Jowelllnq  (Figure  26).  The 
shallowing  prescribed  by  equations  (10)  and  (11)  due  to  the  excessive 
heating  tends*  Initially,  to  override  any  effect  of  the  upwelllng. 

By  hour  six  there  Is  only  a  .4  meter  depth  difference  over  that  of  S2, 
and  essentially  no  temperature  difference  In  the  mixed  layer.  After 
that  point,  the  upwelllng  begins  to  combine  with  the  mixing  to  augment 
the  cooling  of  the  layer,  and  slowing  Its  deepening. 

The  net  result  after  twenty-four  hours  Is  to  produce 
a  layer  that  Is  about  as  deep  as  that  of  S2*  but  that  Is  as  cool  as 
the  case  of  the  evening  upwelllng.  This  represents  a  substantial 
energy  difference  in  the  runs.  Again,  a  slight  surface  temperature 
signature  could  be  seen  with  this  run  In  comparison  with  the  standard. 

3,  Case  III 

The  next  test  Is  that  of  the  Interaction  of  mixing  with  a 
single  cycle  Internal  wave.  By  stepping  the  phase  of  an  Input  sinusoid 
through  the  twenty-four  hour  progression  of  the  model,  a  quasi  two- 
dimensional  motion  can  be  prescribed. 

For  a  wave  of  the  form 

wave  amplitude  cc  A  sin  (wt  ♦  kx) 
the  vertical  motion  would  be  of  the  form 

w(x,z,t)  cc  /V*cos(*t  +  kx)  (21) 

as  seen  In  Figure  6.  By  Increasing  the  value  of  kx,  the  wave  is 
made  to  travel  In  the  -x  direction. 
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The  Initial  conditions  used  were  those  of  cases  SI  and  S2.  The 
vertical  motion  was  Introduced  at  hour  one,  and  the  amplitude  was  varied 
as  t  varied  In  equation  (21),  to  give  a  wave  Induced  vertical  motion. 
This  was  then  run  twenty-four  times,  each  time  starting  at  a  different 
phase  of  the  wave,  to  test  the  Interrelation  of  the  wave  Input  with  the 
cyclical  heat  flux.  For  simplicity,  only  four  major  phase  relation¬ 
ships  are  discussed  here,  as  they  will  demonstrate  the  primary  effects. 
Also,  here,  the  two-dimensional  field  will  be  discussed  for  the  first 
time. 

The  wave  form  of  the  desired  vertical  motion  does  appear  In  all 
four  of  the  displayed  runs.  The  difference  In  each  run  Is  the  phase  of 
the  vertical  motion  in  relation  to  the  heat  flux  (dally  solar  heating 
cycle). 

a.  Shallow  Initial  Layer  (Case  Ilia)  (Figures  27-33) 

%  In  all  four  phase  relationships  presented,  maximum  downward 
motion  at  hours  1,  7,  13,  and  19,  the  Initial  conditions  are  those  of 
SI.  The  magnitude  of  the  peak  upward  and  downward  motions  Is  36  cm/hr 
at  ten  feet,  as  before.  The  wave  Is  Introduced  at  hour  one,  and  pro¬ 
gresses  through  the  system  through  hour  twenty-four. 

(1)  Peak  Upward  at  Hour  7  (Figure  27).  In  this  case  the 
Initial  upward  motion,  hours  one  through  twelve,  the  heating  and  the 
surface  wind  stress  combine  to  give  a  slowly  deepening  layer,  which, 
due  to  the  upward  advectlon  of  the  layer  and  of  cooler,  below  layer 
water.  Is  almost  .1*C  cooler  at  Its  maximum  temperature  than  SI.  At 
the  end  of  the  upward  cycle  the  layer  has  deepened  about  two  meters 
less  than  SI  at  the  same  time,  and  has  set  up  a  slightly  larger  AT  at 
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the  entrainment  zone.  As  a  result*  the  following  downward  cycle  pro¬ 
duces  a  smaller  effect  In  combination  with  the  wind  stress  In  deepen¬ 
ing  than  that  seen  previously.  Because  of  the  Initial  upward  motion 
effects  on  the  temperature  profile,  the  combined  effect  of  wind  stress 
mixing,  surface  cooling,  and  the  downward  motion  produce  a  layer  that, 
after  twenty- four  hours.  Is  only  two  meters  deeper  than  SI.  Addition¬ 
ally,  the  downward  effect  of  the  motion  on  the  temperature  gradient 
Is  sufficient  only  to  make  up  about  60  percent  of  the  deviation  In 
temperatures  caused  by  the  upwelllng  motion,  as  to  that  of  SI. 

The  maximum  #T  for  this  case  Is  .245#C.  It,  also, 
produced  the  deepest  layer  of  the  four  runs  discussed  here,  over  the 
twenty-four  hour  period. 

(2)  Peak  Upward  at  Hour  1  (Figure  28).  In  this  case  the 
upward  motion  was  present  during  the  first  six  hours  of  the  heating 

cycle  only.  The  downward  motion  cycle  then  progressed  for  twelve  hours, 

* 

followed  by  six  hours  of  upward  motion  to  complete  the  twenty-four 
hour  sequence. 

During  the  first  six  hours  the  layer  deepened  to  within 
less  than  one  meter  of  the  depth  reached  by  hour  six  In  SI,  but  the 
layer  temperature  Is  seen  to  be  about  .01"C  cooler.  The  reversal  of 
the  motion  was  too  late  In  the  heating  cycle  to  overcome  that  trend, 
and  by  hour  nine,  the  time  of  the  maximum  surface  temperature  for  Si, 
this  experiment  produced  a  surface  temperature  that  Is  ,013*C  cooler. 

The  layer  depth  reached  here  Is  within  half  a  meter  of  the  SI  depth 
for  the  same  time.  The  difference  Is  attributed  to  the  initial  lift¬ 
ing  of  the  entrainment  zone. 
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As  the  downward  nation  continues*  along  with  the  sur¬ 
face  cooling*  the  layer  deepens  more  quickly  while  the  cooling  of  the 
layer  slows  as  the  effect  of  mixing  Is  reduced.  A  maximum  depth  Is 
reached  two  hours  after  the  downward  motion  has  stopped,  when  the 
upward  motion  Is  of  sufficient  magnitude  to  overcome  the  wind  stress 
mixing.  The  layer  becomes  shallower  from  Its  extreme  depth  of  15.56 
meters  to  14.94  meters  by  hour  twenty-four*  with  associated  layer 
cooling  due  to  surface  heat  loss,  and  Increased  entrainment  due  to 
the  Increased  effect  of  the  wind  stress  mixing. 

The  maximum  ST  for  this  run  Is  .189°C.  The  final 
layer  depth  Is  within  .3  meters  of  SI  for  hour  twenty-four.  The 
variations  of  depth  and  temperature  during  the  run  are  obviously 
much  greater  than  for  SI. 

(3)  Peak  Upward  Motion  at  Hour  19  (Figure  29).  The  fact 
that  the  downwelllng  and  the  heating  cycle  are  In  phase,  as  In  Case 
Ila*  produces,  by  hour  nine,  a  layer  that  Is  two  meters  deeper  and 
.006*C  warmer  than  SI  for  this  same  time.  The  following  upward  motion 
Is  sufficient  to  halt  layer  deepening  at  hour  sixteen  in  opposition 

to  the  mixing,  but  Is  not  sufficient  to  cool  the  layer  as  much  as  for 
SI.  The  final  layer  temperature  Is  .035°C  warmer  than  that  in  SI, 
while  the  final  layer  depth  Is  about  one  meter  less,  though  at  hour 
sixteen  a  depth  of  14.2  meters  was  reached.  The  maximum  ST  for  the 
mixed  layer  In  this  case  was  .181*c. 

(4)  Peak  Upward  Motion  at  Hour  13  (Figure  30).  This 
case  produced  the  warmest  layer  temperature  at  hour  nine,  though  the 
downward  motion  was  only  present  during  the  first  six  hours  of  the 
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surface  heating  cycle.  The  Initial  downward  motion  deepened  the  layer , 
and  thereby  reduced  mixing  In  this  period,  over  the  SI  case,  allowing 
the  surface  layer  to  be  heated  to  14.988'C,  an  Increase  of  .Q1*C  over 
that  In  SI.  The  following  upward  motion  was  sufficient,  hours  eleven 
through  fourteen,  to  hold  the  mixed  layer  at  a  relatively  constant 
depth.  This  allowed  a  relatively  constant,  and  strong  wind  stress 
entrainment  to  combine  with  surface  cooling  to  rapidly  cool  the  layer. 
The  result  after  twenty-four  hours  is  to  produce  a  layer  of  about  the 
same  depth  as  that  of  Si,  but  one  that  is  slightly  cooler. 

The  9T  here  was  .237*0,  a  very  noticeable  difference 
from  SI,  though  the  final  product  did  not  vary  much  In  layer  temperature 
or  layer  depth. 

(5)  Two  Dimensional  Data  (Figure  31,  32  and  33).  The  two 
dimensional  oceanic  field  was  produced  as  previously  discussed.  The 
data  covered  here  Is  that  of  the  evolution  of  the  ocean  system  at 
hours  6,  12,  18,  and  24. 

By  hour  six  the  wave  has  progressed  Into  the  oceanic 
field,  and  Its  effect  below  the  layer  Is  dramatic,  (Figure  31a).  It 
has  produced  a  mixed  layer  that  has  a  nean  depth  of  about  9  meters 
(Initial  h  was  4,5  meters)  with  a  mean  temperature  of  14.957*C.  A 
temperature  structure,  as  could  be  measured  at  the  surface  In  x,  is 
being  formed.  The  present  maximum  difference  In  surface  temperatures 
(In  x)  Is  .016*0. 

By  hour  twelve  (Figure  31b)  the  effect  of  the  wave  Is 
at  Its  maximum,  as  seen  In  the  amplitude  of  the  contours  below  the 
layer.  The  layer  Is  generally  deeper  than  at  hour  six  with  a  mean 
depth  of  about  11.5  meters.  The  deviation  from  the  mean  Is  much 
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greater  for  this  tine,  however.  The  mean  of  the  layer  temperature 
has  gone  to  14.9S0°C,  a  slight  decrease  from  hour  six,  but  more 
significant  Is  the  deviation  of  the  surface  temperatures  In  the  x 
direction.  There  Is  now  a  very  definite  variation  In  x  of  .Q51°C. 

The  layer  has  deepened  to  a  mean  depth  of  about  13 
meters  by  hour  eighteen  (Figure  31c),  with  the  effect  of  the  wave 
beginning  to  be  reduced  in  the  system.  The  mean  layer  temperature 
has  dropped  to  14.846*C,  a  cooling  of  almost  .05*C  in  the  six  hours 
between  Illustrations.  The  surface  temperature  deviation  has  reached 
.091°C.  Warm  and  cold  temperature  bands  are  becoming  defined  In  x. 

The  single  wave  has  passed  through  the  ocean  by  hour 
twenty-four,  and  only  the  remnant  of  Its  passage  Is  seen  (Figure  31a). 
In  an  ocean  that  Is  defined  by  SI  conditions  everywhere,  the  twenty- 
four  hour  product  would  be  an  ocean  with  a  mixed  layer  depth  constant 
at  14.65  meters,  and  mixed  layer  temperature  of  14.767°C.  With  the 
Introduction  of  a  single  diurnal  wave  cycle,  a  variable  condition 
ocean  has  been  produced.  The  mean  depth  of  this  ocean  is  15.012 
meters,  with  local  variations  of  up  to  +1.73  meters  and  -1.32  meters, 
with  depths  at  small  x  being  somewhat  greater  than  depths  at  large  x. 
The  mean  layer  temperature  Is  14.764°C,  but  more  Importantly,  the 
surface  temperature  variation  Is  now  .095°C,  In  a  clearly  defined 
pattern.  The  pattern  that  exists  Is  that  of  a  narrow  zone  of  rela¬ 
tively  cool  water  centered  near  lex  ■  5,  with  a  wider  zone  of  relatively 
warmer  water  centered  near  kx  ■  18,  Below  the  layer  the  Isotherms  are 
found  to  be  parallel  after  the  passage  of  this  single  wave  cycle, 
b.  Deep  Initial  Layer  (Case  Illb)  (Figures  34-40) 

The  Initial  conditions  for  these  runs  were  the  same  as 
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for  S2.  The  wave  cycle  passed  through  the  system  as  In  the  previous 
runs  (Ilia).  The  wave  form  Is  harder  to  discern  In  these  runs,  but 
looking  at  the  contours  below  the  original  layer  (Figures  34-37)  the 
sane  form  as  the  previous  run  are  apparent. 

(1)  Peak  Upward  {lotion  at  Hour  7  (Figure  34).  As  in  the 
previous  deep  initial  layer  cases,  the  shallowing  due  to  the  surface 
heating  drives  the  model  in  the  first  six  hours.  The  coincident  up¬ 
ward  motion  during  this  time  combines  with  the  shallowing  and  mixing 
to  produce  a  layer  that  is  about  1.S  meters  shallower,  but  at  the 
same  temperature  as  In  S2,  by  hour  ten.  Two  noticeable  effects  during 
this  period  of  upward  motion  are  (1)  a  slight  shallowing  of  the  newly 
formed  surface  layer,  and  (11)  the  upward  transport  of  cooler  below 
layer  water,  setting  up  a  larger  AT  at  the  entrainment  zone  for  this 
case  over  that  of  the  standard.  The  downward  motion  then  contributes 
to  the  deepening  of  the  layer  as  does  the  mixing  which  must  initially 
act  over  a  relatively  shallow  depth.  This  entrains  cool  water,  form¬ 
ing  a  quickly  deepening  and  cooling  surface  mixed  layer.  As  the  down¬ 
ward  motion  decreases  the  cooling  slews,  and  the  Increase  in  the  depth 
of  the  layer  Is  seen  to  be  caused  mostly  by  the  vertical  motion.  Of 
all  the  phases  examined  this  produced  the  deepest  layer  (30.34m)  and 
has  the  second  coolest  layer  at  the  end  of  the  period  (15.059*0  in 
relation  to  the  other  three  cases  discussed. 

The  8T  for  this  case  Is  .115°C.  This  Is  about  half 
that  for  the  shallow  layer  case  discussed  previously.  The  original 
deep  Interface  can  be  followed  on  the  contours  (Figure  34b)  over  the 
twenty-four  hour  period.  During  the  downward  motion  It,  again.  Is 
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diffusing,  while  during  the  upward  case  there  Is  a  slight  reinforce¬ 
ment  of  the  Interface. 


(2)  Peak  Upward  Motion  at  Hour  1  (Figure  35).  Again  It 
Is  seen  that  the  shallowing  mode  of  the  mixed  layer  model,  equations 
(10)  and  (11),  controls  the  layer  development  to  hour  ten.  The  layer 
depth  reached  through  the  shallowing,  upward,  and  downward  motion 
combination  by  hour  ten  has  reached  a  depth  that  is  about  .5  meters 
deeper  than  that  In  S2  for  the  same  time.  The  temperature  is  the  same 
as  S2.  Also,  the  downw^/d  motion  makes  entrainment  of  water  Into  the 
mixed  layer  due  to  mixing  more  difficult,  such  that  cooling  after  hour 
ten  progresses  much  more  slowly  than  In  that  of  the  previous  case,  and 
that  of  S2.  The  net  result  In  the  layer  depth  of  this  run  Is  a  layer 
that  Is  about  as  deep  as  that  of  S2  after  the  twenty-four  hour  period. 
The  model  did  not  deepen  more  than  S2  due  to  the  dominant  upward  motion 
during  the  last  six  hours.  The  mixed  layer  temperature  was  slightly 
warmer  than  that  of  S2  at  the  end  of  the  run  (.071°C).  Also,  it  can 
been  seen  that  the  below  layer  water  is  somewhat  warmer  than  In  S2  for 
the  end  of  the  period,  due  to  the  combined  effects  of  the  heating  and 
downward  vertical  motion  timing. 

The  ST  for  this  case  Is  .092°C.  This  is  less  signi¬ 
ficant  than  In  the  previous  run. 

(3)  Peak  Upward  Motion  at  Hour  19  (Figure  36).  Again, 
the  shallowing  mode  of  the  model  controls  the  Initial  ten  hour  layer 
temperature  adjustment.  The  downward  vertical  motion,  however,  drives 
the  layer  two  meters  deeper  during  this  time  causing  reduced  mixing 
due  to  the  surface  wind  stress.  The  upward  motion  keeps  the  layer 
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from  deepening  as  quickly  as  before  after  hour  twelve.  Because  of 
the  below  layer  temperature  adjustment  by  the  downward  motion  and 
the  Initial  deepening  of  the  layer,  the  layer,  also,  does  not  cool 
as  quickly.  This  case  produced  the  shallowest  final  layer  depth, 
and  the  second  wannest  layer  at  the  end  of  the  period,  relative  to 
the  other  three  runs. 

The  8T  for  this  case  Is  .103°C,  which  Is  about  the 
same  as  the  first  case  which  was  180°  out  of  phase  with  this  run. 

(4)  Peak  Upward  Motion  at  Hour  13  (Figure  37) .  The 
coincident  shallowing,  wind  stress  mixing,  heat  flux,  and  cycled 
vertical  motion  over  the  first  ten  hours  produce  a  surface  mixed 
layer  that  has  the  same  temperature  as  S2,  and  almost  the  same  layer 
depth  (to  within  less  than  .5  meters).  Below  the  layer,  however, 
the  upward  motion  has  made  relatively  cooler  water  available  by  hour 
ten  (In  comparison  with  S2),  and  the  layer,  though  not  rapidly  deep- 
enlng  due  to  the  continued  upward  motion,  cools  quickly  as  large 
quantities  of  water  are  entrained  over  a  relatively  shallow  depth, 
a  trend  which  continues  to  the  end  of  the  period.  The  net  result  Is 
a  mixed  layer  that  Is  1.5  meters  deeper  than  that  of  S2,  but  Is  also 
the  coolest  layer  of  the  four  runs.  Because  of  the  below  layer  temp¬ 
erature  profiles  which  are  set  up  by  the  upward  motion  for  the  mid¬ 
twelve  hours,  the  continued  run  of  this  model  over  a  second  day  gives 
continued  rapid  cooling  of  the  layer  due  to  Increased  entrainment  of 
cooler,  upward  advected  water,  and  to  Increased  effects  of  mixing. 

A  noticeable  deepening  does  take  place  at  the  end  of  the  period  due 
to  the  downward  motion,  which  accounts  for  the  fact  that  this  run 
ended  deeper  than  S2. 


67 


The  ST  for  this  run  is  ,125°C.  This  is  the  largest 
of  the  four  runs,  but  it  Is  still  only  half  of  the  maximum  8T  seen  for 
the  shallow  initial  layer  case. 

(5)  Two-Dimensional  Data  (Figure  38.  39  and  40).  The  two- 
dimensional  oceanic  field  Is  produced,  here,  as  previously  discussed, 
for  the  single  wave  cycle.  The  data  covered  is  for  hours  6,  12,  18, 
and  24.  A  common  feature  is  the  varying  below  layer  temperature  struc¬ 
ture  set  up  because  of  the  Introduction  of  the  initial  deep  layer. 
Because  of  this,  the  changes  in  the  thermocllne  below  the  initial  deep 
Interface  best  show  the  effect  of  the  wave  progression  through  the 
system. 

As  In  the  shallow  case,  by  hour  six,  the  wave  has 
progressed  into  the  oceanic  system  producing  some  horizontal  sea  sur¬ 
face  temperature  structure.  However,  the  effect  on  the  surface 
temperature  Is  almost  negligible.  In  that  the  variability  of  the 
surface  temperature  In  x  Is  only  ,001"C.  This  Is  due  to  the  shallow¬ 
ing  mode  response  of  the  model  over  the  first  six  hours.  Even  with 
this,  however,  the  depth  structure  has  varied  (Figure  38a)  so  that 
there  has  been  a  change  from  the  constant  structure  that  would  be 
seen  In  $2  of  11.767  meters  for  the  layer  depth,  to  a  layer  of  mean 
depth  of  11.769  meters  t  >35  meters,  indicating  that  an  adjustment 
due  to  the  wave  Is  taking  place. 

By  hour  twelve  much  more  temperature  and  layer  (h(x)) 
variation  Is  beginning  to  take  shape.  At  hour  twelve  the  case  without 
the  vertical  motion  has  a  layer  depth  of  15.037  meters,  and  a  layer 
temperature  of  15.159°C  everywhere.  In  this  case,  where  the  maximum 
effect  of  the  wave  has  been  reached,  the  layer  depth  has  a  mean  of 
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15.231  meters  with  a  deviation  of  t  2.5  meters.  The  variability  In 
the  surface  temperature  In  x  has  reached  .01*C.  Though  the  variabil¬ 
ity  In  the  surface  temperature  Is  not  very  large,  a  pattern  Is  again 
beginning  to  emerge  In  the  x  direction.  At  small  x  cooler  Mater  Is 
forming,  Mhlle  at  large  x  a  band  of  warmer  water  Is  emerging. 

The  layer  has  continued  to  deepen  to  hour  eighteen. 
The  layer  depth  has  Increased  to  a  mean  of  19.46  m,  t  3.4  m,  and  the 
variation  In  the  surface  temperature  structure  In  x  has  reached  .03*C. 
The  warm-cold  pattern  seen  at  hour  six  Is  repeated  and  strengthened. 
Though  the  temperature  difference  between  warm  and  cool  water  Is  not 
as  great  as  In  the  shallow  case,  the  banding  of  the  patches  Is  much 
more  distinct  here  (Figure  38,  Figure  31). 

At  the  end  of  the  period  the  wave  has,  again,  been 
removed  from  the  system,  and  only  the  remnant  of  Its  passage  Is  seen. 
The  final  layer  depth  Is  deeper  than  In  the  previous  hours,  and  much 
deeper  than  In  the  shallow  case.  The  variability  of  h  In  x  Is  also 
much  greater  here  than  In  the  shallow  case.  The  surface  temperature 
patchiness  Is  well-defined,  but  the  variability  of  the  temperature  In 
x  Is  only  ,035°C. 

Comparing  the  two  cases,  shallow  and  deep  Initial 
layers.  It  Is  seen  that  for  the  same  amplitude  wave  the  model  produces 
two  distinct  ocean  types.  For  the  shallow  case  an  ocean  Is  formed 
that  has  a  relatively  shallow  layer  with  a  large  temperature  vaca¬ 
tion  In  the  x  direction,  but  with  that  temperature  variation  primarily 
in  wide  bands.  For  the  deep  case  a  deep  surface  layer  Is  formed  with 
smaller  variability  In  the  surface  temperature.  However,  the  temp-, 
erature  structure  In  x  Is  much  better  defined.  Also  In  the  deep  case. 
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the  patches  of  warm  and  cold  water  are  much  closer  to  each  other  (In  x). 
This  could  suggest  a  mechanism  for  generating  or  sharpening  SST  fronts 
through  surface  heating,  and  Internal  solitary  wave  passage. 

4.  Case  IV 

It  was  noted  In  Case  III  that  after  twenty-four  hours  the  wave 
form  had  passed  through  the  ocean  system,  and  had  left  behind  only  the 
residue  of  Its  effect  on  the  system,  over  and  above  the  standard  effects. 
This  can  be  compared  to  a  solitary  wave  passing  through  an  ocean  system. 
While  this  has  theoretical  advantages.  It  Is  not  a  very  realistic  prob¬ 
lem,  In  that  wave  forms  seldom  take  on  this  sort  of  motion.  One  pos¬ 
sible  exception  may  be  solltons  generated  at  extreme  tidal  stages 
(Apel,  1982). 

The  next  progression,  therefore  Is  to  look  at  a  continuous  wave 
system  that  passes  through  the  ocean  with  constant  form,  until  the  run 
Is  terminated.  This  Is  done  by  simply  modifying  the  mathematics  of 
the  previous  case  to  Include  initialization  of  the  temperature  and 
salinity  profiles  with  the  wave  form  according  to  the  following. 

We  know  that 

vKx,x,t)  m  A^z.co»(wt  +  kx) 
and 


8b  m  ~f  wdl 

therefore 

8k  «  -  Ayiiwtwt  ♦  kx) 
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For  initial  displacement 


8h(x,z,0)  m  fth0 


therefore 


tll0  ■ 


-A^alw(kx) 


(22) 


Applying  equation  (22)  to  form  the  initialized  layer  depth,  and  com¬ 
bining  this  with  equation  (16)  to  initialize  the  temperature  profile, 
a  continuous  wave  field  can  now  be  passed  through  the  model  ocean 
(Figure  7). 

From  this  point  on  the  discussion  will  be  confined  to  the  deep 
Initial  layer  case,  S2  standard,  and  variations  thereof.  This  is  done 
because  all  the  primary  effects  discussed  to  this  point  appear  In  both 
the  shallow  and  deep  cases,  as  do  the  effects  to  be  discussed  in  vary¬ 
ing  degrees,  and  It  Is  not  necessary  to  continuously  reiterate  this 
discussion. 

a.  Diurnal  Wave  (Case  IVa)  (Figures  41,  42  and  43) 

Because  of  the  overriding  effect  of  the  shallowing  mode  of 
the  model  In  the  first  six  hours  of  the  runs,  the  results  produced 
here  were  exactly  the  same  as  those  discussed  In  Case  Illb,  for  the 
phase  relationships.  This  was  of  course  expected  from  the  theory. 

(1)  Two-Dimensional  Data  (Figure  44.  4S  and  46).  Again, 
because  of  the  dominance  of  the  shallowing  mode,  the  hour  6,  12,  18, 
and  24  two-dimensional  fields  produced  here  are  the  same  as  those  of  * 
Case  Illb  for  the  area  of  the  entrainment  zone,  and  the  mixed  layer. 
Below  the  layer,  however,  the  wave  form  is  now  seen  to  progress  through 
the  system  In  a  continuous  mode.  This  is  best  seen  In  the  contours 
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Figure  7.  Diurnal  wave  relationship  (of  continuous  diurnal  waveform  to  model  location  (x) 
Amplitude  variable.) 


In  which  the  amplitude  of  the  system  remains  a  constant,  but  the  wave 
form  Is  shifted  with  time. 

% 

It  can  be  seen  that  the  patch  of  warmer  waters  forms 
(Figure  44)  and  takes  a  phase  that  barely  leads  the  peak  of  the  wave 
system.  The  entrainment  zone  (  h(x)  )  now  takes  the  form  of  the 
passing  wave  system,  as  would  be  expected,  and  yet  was  not  obvious  In 
the  previous  case.  The  divergence  and  convergence  of  below  layer 
fields  Is  also  obvious  here,  where  previously  It  was  not  particularly 
obvious. 

b.  Semi-diurnal  Wave  (Case  IVb)  (Figures  47-50) 

A  dominant  Internal  wave  system  of  the  temperate  oceans 
Is  the  semi-diurnal  tidal  cycle.  South  of  30°  N  there  Is  also  a 
strong  diurnal  tidal  cycle.  A  semi-diurnal  tidal  field  can  be  pro¬ 
duced  for  testing  In  this  model  by  simply  adjusting  the  frequency  of 
the  prescribed  w  for  equation  (21),  and  setting  the  Initial  profiles 
with  the  appropriate  frequency  of  equation  (22)  and  (17). 

The  amplitude  of  the  motion  remains  the  same  as  the 
diurnal  cycle  case.  Because  the  cycle  of  wave  mixing  interactions 
repeats  Itself  after  a  2r  phase  shift,  which  occurs  after  12  hours 
In  this  case,  only  the  cases  of  six  and  twelve  hour  phase  shifts 
will  be  discussed.  The  18  and  24  hour  phase  shifts  are  merely  repeats. 
(1)  Six  Hour  Shift  (Figure  47a  and  b).  Figure  8  shows 

e 

the  relationship  of  the  heating  cycle  to  the  wave  motion  In  this  case. 
The  vertical  motion  Is  seen  to  Imaedlately  combine  with  the  shallowing 
mode  of  the  model  through  the  Initialization  to  give  a  somewhat 
shallower  Initial  layer.  The  shallowing  mode  of  the  model  holds  the 
layer  temperature  the  same  as  In  S2  to  hour  six,  while  the  vertical 
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Figure  8.  Phase  relationship  of  Case  IVb  (six  hour  shift)  of  semi-diurnal  waveform  to  diurnal 
surface  heating  cycle  over  time  (hrs).  Amplitudes  are  not  proportional. 


motion  varies  the  layer  depth,  and  alters  the  below  layer  temperature 
structure.  By  hour  six  a  half-cycle  of  the  wave  has  passed  through 
the  model  point,  and  a  layer  has  been  produced  that  Is  almost  one 
meter  deeper  than  S2  at  this  time,  but  which  has  the  same  layer  temp¬ 
erature.  Continued  but  decreasing  downward  motion  combines  with  the 
downward  heat  flux  until  hour  nine  to  slowly  deepen  the  layer,  but 
the  layer  continues  to  warm.  This  process  continues  through  hour  ten, 
though  upward  motion  has  begun  because  the  heating  and  entrainment  of 
relatively  warm  water  are  greater  than  the  effect  of  the  resistance  to 
deepening  due  to  the  surface  heat  flux,  and  the  small  upward  motion. 

By  hour  eleven  die  upward  motion  dominates,  and  the  layer  shallows 
while  It  cools  slightly  due  to  the  Increased  entrainment.  At  hour 
twelve  the  layer  Is  as  shallow  as  It  will  ever  get,  and  the  first  wave 
has  passed,  flow  decreased  upward  motion  followed  by  downward  motion, 
combined  with  surface  cooling  and  mixing  through  a  progressively  cooler 
surface  mixed  layer  quickly  Increases  the  layer  depth  to  a  maximum  of 
26.7  meters  at  hour  twenty-one.  This  Is  deeper  than  S2  gets  In  the 
entire  twenty-four  hour  period.  Yet,  the  layer  Is  not  as  cool  as  In 
the  S2  case.  By  hour  twenty-one  the  upward  motion  Is  again  dominant 
and  the  layer  shallows. 

(2)  Twelve  Hour  Shift  (Figure  47c  and  d).  Figure  9  shows 
the  relationship  of  the  heating  cycle  with  the  wave  In  this  case.  The 
Initialization  combines  with  the  vertical  motion  and  shallowing  mode 
to  produce  an  initial  layer  that  Is  deeper  than  that  In  S2,  but  soon 
becomes  shallower,  as  the  new  surface  layer  Is  established.  By  hour 
six  a  layer  of  the  same  temperature  as  S2,  but  .72  meters  shallower 


has  been  formed.  The  net  result  of  the  cyclical  interactions  In  this 
case  Is  a  layer  that  Is  1.5  meters  deeper  than  In  S2,  but  which  is  at 
about  the  same  temperature.  The  difference  In  the  layer  temperature 
between  this  run,  and  the  previous  Is  .004*C.  However,  the  layer 
depth  Is  2.2  m  different. 

(3)  Two-Dimensional  data  (Figure  48.  49  and  501.  At  hour 
six  a  very  dlscemable  depth  structure  has  been  set  up  at  the  bottom 
of  the  mixed  layer  which  reflects  the  structure  of  the  passing  wave. 

Above  the  entrainment  zone  the  temperature  structure,  as  produced  by 
the  shallowing  mode,  is  homogeneous.  Below  the  layer  the  structure  Is 
highly  varied  In  x  reflecting  the  combined  effects  of  the  vertical 
motion  to  this  time.  The  mean  layer  depth  Is  11.75  meters,  while  that 
of  S2  Is  11.77  m,  showing  the  domination  of  the  shallowing  to  this  point. 
A  variation  of  h  In  x  of  *  .75  meters  Is  already  exhibited. 

By  hour  twelve  the  mixed  layer  depth  structure  In  x, 
and  the  layer  temperature  variability  In  x  have  evolved  considerably 
over  that  of  hour  six.  At  the  entrainment  zone  the  structure  Is  approx¬ 
imately  sinusoidal,  as  expected.  Over  the  last  six  hours  a  large 
amount  of  water  from  below  the  layer  has  been  entrained  while  surface 
heating  has  occurred,  and  the  mixed  layer  has  heated  and  cooled 
accordingly  to  begin  to  set  up  the  pattern  of  warm  and  cold  bands  of 
temperature  in  x.  The  mean  depth  of  the  layer  Is  now  15.21  meters 
t  2.50  meters,  while  that  for  S2  Is  15.04  maters.  The  temperature 
variation  In  x  as  seen  at  the  surface  Is  now  .012*C. 

The  structure  of  tempera  tuna  and  layer  depth  profile 
has  been  further  accented  by  hour  eighteen.  The  layer  has  continued 
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to  deepen,  and  the  fluctuation  in  x  Increased.  The  banding  of  the 
In  layer  temperature  pattern  In  x  has  become  more  defined. 

After  twenty* four  hours  a  layer  has  been  produced 
that  has  a  depth  of  26.14  m,  t  3.3  m.  The  surface  temperature  pattern 
In  x  Is  quite  clearly  defined,  though  the  variation  in  temperature  Is 
only  .015*0. 

5.  Case  V  (Superimposed  Waves)  (Figure  51.  52  and  S3) 

Since  oceanic  Internal  waves  are  almost  never  simple  sinusoids, 
one  more  progression  Is  considered.  Spectral  analysis  of  the  MILE 
data  (Figure  10}  at  station  P  (55*  N)  shows  the  dominant  wave  field  to 
be  the  semi-diurnal  tide  superimposed  by,  at  least,  the  6*7  hour  har¬ 
monic  of  the  tidal  wave.  Internal  wave  theory  Indicates  that  the 
period  of  the  Internal  waves  must  be  less  than  or  equal  to  the  Inertial 
period.  North  of  approximately  30*  N,  no  diurnal  internal  waves  may 
propagate,  and  therefore,  they  are  not  seen  In  the  MILE  data.  South 
of  this,  however,  diurnal  and  semi-diurnal  tidal  period  waves  may 
propagate.  Therefore,  since  previous  discussions  have  been  concerned 
with  diurnal  and  semi-diurnal  waves,  and  there  Is  justification  for 
their  presence,  these  waves  will  be  superimposed  to  form  a  slightly 
mere  complicated  wave  system. 

The  combined  wave  Is  directed  through  the  two-dimensional 
ocean  which  has  the  Initial  conditions  of  S2,  from  right  to  left 
(large  kx  to  small  kx).  The  amplitude  of  the  motion  Is  such  that 
both  waves  peak  at  the  same  amplitude  at  their  peak  motions,  though 
one  has  half  the  period  of  the  other.  Both  waves  start  with  an 
amplitude  of  vertical  motion  In  phase  at  time  one. 
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The  superimposed  waves  reflect  all  of  the  Interactions  discussed 
In  the  previous  sections.  Further,  this  simple  addition  of  two  wave¬ 
forms  has  developed  several  new  features.  The  new  features  can  be 
traced  In  their  development  from  hour  one  through  hour  twenty- four, 
where  they  become  most  evident.  Since  this  Is  the  case,  the  discussion 
will  be  limited  to  the  case  of  the  two-dimensional  ocean  at  hour  twenty- 
four. 

By  this  time  a  complex  structure  has  emerged  for  the  two-dimen¬ 
sional  ocean,  both  In  temperature  and  In  profile  of  the  layer  depth. 

An  ocean  has  been  formed  that  has  a  mean  layer  depth  of  26.78  ra,  but 
this  depth  has  a  variation  that  Is  not  symmetrical  from  small  to  large  x. 
At  small  x  the  variation  from  the  mean  depth  Is  seen  to  be  a  maximum  of 
-5.45  m,  while  at  large  x  the  layer  depth  varies  to  a  depth  deviation 
from  the  mean  of  +7.83  m.  An  ocean  has  been  formed  that  has  a  deep 
layer  at  one  end,  and  a  shallow  layer  at  the  other.  Further,  a  devia¬ 
tion  of  the  surface  temperature  In  the  x  direction  of  .042°C  is  seen. 

The  pattern  of  the  w^rm  and  cold  bands  again  Is  seen,  but  this  also  Is 
no  longer  symmetrical.  The  temperature  pattern,  as  seen  at  the  surface. 
Is  a  cold  band  beginning  at  small  kx,  progressing  through  a  weak  gradi¬ 
ent  to  a  warm  band  centered  near  kx  »  15.  A  strong  gradient  Is  then 
encountered  leading  to  a  cold  band  centered  about  kx  »  21.  Through  the 
introduction  of  a  simple  two  wave  system,  an  extremely  diverse  ocean 
has  been  produced. 

Shifting  the  phase  relationships  of  the  two  Internal  waves,  and 
running  the  model  again,  resulted  In  similar  patterns.  Positions  of 
the  temperature  maxima  and  minima  In  the  layer,  and  the  position  of 
the  depth  variations  shift,  but  their  relative  positions  are  maintained. 
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The  temperature  variation  In  the  x  direction  as  seen  at  the  surface 
does  not  change  with  the  phase  changes. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  analysis  of  the  mixing  response  to  the  prescribed  vertical 
motion  shows  that  vertical  advectlon  can  be  successfully  Introduced 
into  this  one* dimensional  mixed  layer  model  without  numerical  diffi¬ 
culties.  This  Is  Important  In  that  It  Indicates  that  a  single  station 
forecasting  model  can  be  altered  to  Include  the  variability  of  long 
Internal  wave  motion.  This  step  is  the  first  of  a  progressive  devel¬ 
opment  of  a  single  station  forecasting  model  that  may  In  the  future 
Include  other  local  advectlve  effects. 

Using  an  ensemble  of  one-dimensional  stations,  a  two-dimensional 
mixing  response  to  prescribed  internal  wave  motion  was  simulated. 
Although  the  Initial  temperature  field  (without  wave)  and  the  surface 
boundary  conditions  were  both  horizontally  homogeneous,  a  significant 
wave-induced  horizontal  structure  was  created  in  the  T(x,z)  field,  doth 
linear  and  non-linear  effects  have  been  demonstrated  in  the  simple 
single  wave  and  superimposed  wave  cases.  Horizontal  variability  has 
been  shown,  both  In  mixed  layer  depth  and  In  mixed  layer  temperature 
structure.  The  mixed  layer  temperature  develops  a  corrollary  structure 
of  alternating  warm  and  cold  patches.  In  the  case  of  the  superimposed 
waves  of  different  frequencies,  this  patchiness  pattern  can  be  complex, 
and  can  give  rise  to  horizontal  gradients  which  could  be  frontogenetlc. 

There  Is  an  Important  relationship  between  the  phase  of  Internal 
motion  and  the  cyclical  diurnal  heat  flux  boundary  condition.  This  Is 
clearly  demonstrated  In  the  numerical  experiments  In  which  the  phase 
of  the  input  pulse  or  wave  system  has  been  varied. 
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In  the  case  of  superimposed  diurnal  and  semi-diurnal  waves,  the 
horizontal  variation  In  temperature  In  the  mixed  layer  was  about  .04*C 
after  only  a  single  diurnal  cycle.  In  the  same  period,  the  mixed  layer 
depth  varied  by  about  13  m.  The  wave  dispersion  relation  for  a  two  layer 
ocean  at  25°  N  gives  a  horizontal  wave  length  scale  of  from  10  to  100  Km, 
depending  upon  stratification.  This  wave-mixing  Interaction  may,  there¬ 
fore,  be  an  Important  mechanism  for  the  generation  of  sea  surface 
patchiness  as  found  by  Butler  (1981). 

Further  study  must  be  done  In  this  area.  Clearly,  more  complicated 
waveforms  will  produce  more  complex  responses.  The  model  should  be 
tested  with  salinity  changes.  Including  evaporation  and  precipitation. 
Tests  should  be  conducted  on  data  sets  such  as  those  of  MILE.  Compar¬ 
ison  of  the  mixed  layer  model  results  using  MILE  boundary  conditions, 
with  and  without  vertical  advectlon,  as  observed  In  spectrum  analysis, 
should  further  demonstrate  the  importance  of  addition  of  vertical 
motion  to  mixed  layer  models. 

Finally,  the  wave  systems  should  be  tested  In  an  Interactive  mode 
with  the  mixing  of  the  model.  Energy  exchanges,  as  prescribed  by  wave 
theory  and  the  studies  of  Linden  (1975),  and  Kantha  et  al.  (1977) 
should  be  allowed,  and  the  relative  significance  of  the  different 
mechanism  should  be  evaluated.  Internal  vertical  motion  Is  an  Impor¬ 
tant  dynamic  variable,  and  should  be  treated  as  such  In  the  future  of 
mixed  layer  modelling. 
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0.50 
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l 

IS 
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0.54 
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-.0354* 
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l 

14 

11. 5 

14.075 

14.000 

0.50 
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-.01110 

-14.504 

1 

IT 

11.5 

14.050 

34.000 

0.50 

100.0 

-.01477 

-40.742 

L 

14 

11.  S 

14.457 

24.000 

O.so 

100.0 

-. 01400 

-42.100 

1 

1* 

11.4 

14.020 

14.  MO 

0.50 

100.0 

-.01477 

-40.444 

1 

20 

11.7 

14.004 

24.000 

0.50 

100.0 

-.03114 

-34.911 

t 

21 

11.4 

14.707 

14. MO 

o.so 

100.0 

-.02344 

-29.944 

x 

1 

22 

11.7 

14.770 

14. MO 

o.so 

100.0 

-.01400 

-31.220 

K 

l 

21 

11.4 

14.734 

14.  MO 

0.20 

100.0 

-.00412 

-11.10ft 

a 

1 

24 

12.1 

14.737 

34.000 

o.so 

100.0 

0.0 

0.0 

S 

2 

1 

12.4 

14.720 

14.000 

o.so 

-  1.5 

0.0 

o.o 

2 

2 

12.2 

14.727 
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o.so 

-100.0 

0.0 

0.0 

2 

2 

12.7 

14.711 

14.000 

o.so 
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0.0 

0.0 

2 

4 

12.4 

14.741 

14.000 

o.so 
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0.0 

o.o 

2 

S 

12.9 

14.751 

14.000 

0.10 

-100.7 

0.0 

o.o 

2 

4 

11.0 

14.745 

14.000 

O.SO 

-102.5 

0.0 

0.0 

2 

7 

12.1 

14.777 

34.000 

o.so 

-200.4 
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• 

11.2 

14.704 

14. MO 

o.so 

-240.7 

0.0 

o.o 
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9 
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14.7*1 
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0.50 

•1*4.7 

0.0 

0.0 

2 

10 

12.4 

14.7*0 

34.000 

o.so 

-101.2 

0.0 

0.0 

2 
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12.4 

14.704 
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0.0 
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2 

12 

13.4 

14.740 
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0.50 

100.2 

0.0 

o.o 

2 
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13.4 

14.755 
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O.SO 
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2 
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13.4 

14.742 
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0.50 
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2 
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13.4 
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14.000 

O.SO 
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-.02544 
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2 

14 

12.4 
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0.50 

100.0 
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-42.1*9 
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IT 

13.1 
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0.50 
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12.4 
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0.50 
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12.4 
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0.50 

100.0 
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0.0 
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-240.4 
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-104.7 
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0.02940 

29.170 

14.940 
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o.so 

100.2 
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o.so 
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o.so 

100.0 
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o.so 
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20.420 
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0.50 
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14.021 

>4.000 

0.90 
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0.0 

0.0 

14.022 
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o.so 
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0.90 
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14.709 

14.000 

0.90 
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-.01110 
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14.774 
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o.so 
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-.01477 
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o.so 

-  1.9 
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0.90 
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-104.1 
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0.90 

-200.7 

-.01000 

-24.009 

14.010 

>4.000 

O.SO 

•202.9 

-.00912 

-12.003 

14.022 

24.000 

O.SO 

-200.9 

0.0 

0.0 

14.011 

24.000 

0.90 

-240.7 
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0.02477 

00.092 

14.700 

14.000 
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•  l.S 

0.01400 
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*4.41* 

9.7 

14.010 
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l.< 

14.044 

>4.004 

o.so 
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14.0*4 
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14.040 
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19.  4 
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0.0 

0.0 
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-10.04# 

10.0 
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u.ooo 

o.so 

-104.7 

-.01000 

-19.949 

10. • 

14.009 

u.ooo 

o.so 

-191.1 

-.0*944 

•27.9*4 

10.9 

14.97* 

u.ooo 

o.so 

-  4.1 

-.01114 

-14.249 

10.0 

14.090 

u.ooo 

o.so 
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-19.194 
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14.0*7 

*4.000 

o.so 

100.0 
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-*9.42* 

10.7 

14.019 

>4.000 

o.so 

100.0 

-.01477 

-*7.999 

10.7 

14.009 

u.ooo 

0.90 

100.0 

-.01119 

-14.041 

10.0 

14.992 

u.ooo 

o.so 

100.0 

-.oafoc 

-27.090 

u.o 

14.091 

u.ooo 

o.so 

100.0 

-.01900 

-14.993 

11.2 

14.944 

u.ooo 

o.so 
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-.009*1 

-14.499 
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14.014 

u.ooo 

o.so 

100.0 

0.0 

0.0 

u.o 

14.910 

u.ooo 

o.so 

100.0 
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11.119 

u.* 

14. 704 

14.000 

o.so 

100.0 

0.01900 

*2.114 

11.1 

14.770 
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o.so 

100.0 

0.02944 

11.0*9 

ll.t 

14.74* 

u.ooo 

0.10 

100.0 

0.01119 

42.S42 

u.o 

14.7*1 

u.ooo 

o.so 

100.9 

0.02477 
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If.  4 

14. 744 

14.000 

o.so 

-  l.S 

0.01400 

94.497 

U.l 

14.744 

u.ooo 

0.10 

-101.0 

0.01477 
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u.o 

14.7*2 

34.000 

vso 

-104.1 

0.01114 
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17.1 

14.742 

u.ooo 

0  10 

-240.4 

0.0*949 

42.40# 

17.9 

14.774 

u.ooo 

o.so 

-200.7 

0.01900 

11.424 

17.0 

14.704 

u.ooo 

o.so 

-102. S 

0.009*2 

14.101 

17.0 

14.707 

u.ooo 

o.so 

-200.9 

0-0 

0.0 

17.7 

14.007 

u.ooo 

o.so 

-240.7 

-.00932 

-14. SOS 

17.  * 

14.919 

u.ooo 

o.so 

-104.7 

-.01900 

-11.717 

17.1 

14.414 

u.ooo 

o.so 

-102.3 

-.02949 

•44.194 

U.O 

14.410 

u.ooo 

o.so 

-  4.1 

-.0*114 
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U.« 

14.700 

14.000 

o.so 
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-.0*477 

-94.090 

U.O 

14.700 
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o.so 

100.0 

-.0*400 
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10.7 

14.770 

u.ooo 

0.90 

100.0 
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1*.* 

14.740 

14.000 

o.so 

100.0 

-.fiUf 

-49.000 

u.l 

14.797 

u.ooo 

o.so 
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-.03944 

•19.141 

if.  i 

14.740 

u.ooo 

0.99 

100.9 

-.01900 

-27.9*9 
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14.7*1 

14.000 

o.so 

100.9 

-.00912 

-14.472 

11.7 
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0.19 

100.9 

0.0 

o.o 

U.2 

14.700 

14.000 

o.lo 
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0.  099 SI 
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10.7 

14.000 

14.000 

o.so 

100.9 
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17.1 

14.000 

u.ooo 

o.so 

100.0 
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10.0 

14.070 

u.ooo 

o.so 
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